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ABSTRACT 

We  conducted  a  survey  of  nearby  binary  systems  composed  of  main  sequence  stars  of  spectral  types  F  and  G 
in  order  to  improve  our  understanding  of  the  hierarchical  nature  of  multiple  star  systems.  Using  Robo-AO,  the 
first  robotic  adaptive  optics  instrument,  we  collected  high  angular  resolution  images  with  deep  and  well-defined 
detection  limits  in  the  Sloan  Digital  Sky  Survey  V  band.  A  total  of  695  components  belonging  to  595  systems  were 
observed.  We  prioritized  observations  of  faint  secondary  components  with  separations  over  10''  to  quantify  the  still 
poorly  constrained  frequency  of  their  subsystems.  Of  the  214  secondaries  observed,  39  contain  such  subsystems;  19 
of  those  were  discovered  with  Robo-AO.  The  selection-corrected  frequency  of  secondary  subsystems  with  periods 
from  10^-^  to  10^  days  is  0.12  zb  0.03,  the  same  as  the  frequency  of  such  companions  to  the  primary.  Half  of  the 
secondary  pairs  belong  to  quadruple  systems  where  the  primary  is  also  a  close  pair,  showing  that  the  presence  of 
subsystems  in  both  components  of  the  outer  binary  is  correlated.  The  relatively  large  abundance  of  2+2  quadruple 
systems  is  a  new  finding,  and  will  require  more  exploration  of  the  formation  mechanism  of  multiple  star  systems. 
We  also  targeted  close  binaries  with  periods  less  than  100  yr,  searching  for  their  distant  tertiary  components,  and 
discovered  17  certain  and  2  potential  new  triples.  In  a  subsample  of  241  close  binaries,  71  have  additional  outer 
companions.  The  overall  frequency  of  tertiary  components  is  not  enhanced,  compared  to  all  (non-binary)  targets, 
but  in  the  range  of  outer  periods  from  10^  to  10^-^  days  (separations  on  the  order  of  500  AU),  the  frequency  of 
tertiary  components  is  0.16  zb  0.03,  exceeding  the  frequency  of  similar  systems  among  all  targets  (0.09)  by  almost 
a  factor  of  two.  Measurements  of  binary  stars  with  Robo-AO  allowed  us  to  compute  first  orbits  for  9  pairs  and  to 
improve  orbits  of  another  1 1  pairs. 

Key  words:  binaries:  close  -  binaries:  general  -  instrumentation:  adaptive  optics  -  stars:  formation  - 
techniques:  high  angular  resolution 
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1.  INTRODUCTION 

Statistics  of  hierarchical  stellar  systems  provide  important 
clues  to  star  formation  mechanisms  that  are  still  actively 
researched  and  debated.  Current  theories  give  contradictory  pre¬ 
dictions  about  hierarchical  multiples  (for  example,  Y-body  dy¬ 
namics  versus  hydrodynamical  simulations  produce  different 
multiplicity  fractions);  so  far,  none  of  them  is  capable  of  mod¬ 
eling  the  distributions  of  periods,  mass  ratios,  and  hierarchies 
seen  in  observations. 

As  periods  and  separations  of  stellar  pairs  vary  by  several 
orders  of  magnitude,  comprehensive  coverage  of  the  parameter 
space  can  be  achieved  only  by  combining  complementary 
observing  techniques,  and  only  for  the  nearest  stars.  However, 
the  number  of  objects  in  the  best-studied  25  pc  sample  of 
Raghavan  et  al.  (2010)  is  too  small  for  deriving  meaningful 
statistics  of  triple  and  higher-order  multiples.  For  this  purpose, 
a  larger,  distance-limited  sample  of  ~5000  F  and  G  dwarfs 
within  67  pc  of  the  Sun  was  constructed  from  the  Hipparcos 
catalog;  it  is  called  FG-67  hereafter. 

The  sample  definition  and  the  database  are  presented  in 
Tokovinin  (2014).  The  census  of  companions  to  primary  stars  in 
FG-67  is  fairly  complete  over  the  whole  range  of  periods,  with 


an  overall  detection  probability  of  about  80%,  except  for  low- 
mass  companions  at  separations  from  O'.'!  to  20"  (see  details 
in  the  above  paper).  Such  companions  can  be  discovered  by 
high-contrast  imaging  with  adaptive  optics  (AO).  Secondary 
components  can  be  close  pairs  as  well,  but  discovery  of 
subsystems  in  the  secondary  components  is  more  problematic, 
with  average  detection  probability  (before  this  survey)  of  only 
about  15%;  most  subsystems  in  the  secondaries  are  presently 
missed.  If  they  are  as  frequent  as  subsystems  in  the  primaries, 
the  actual  number  of  stellar  hierarchies  in  the  FG-67  sample  is 
substantially  larger  than  known  today. 

It  is  important  to  test  the  frequency  of  subsystems  in  the  secon¬ 
daries  for  constraining  mechanisms  of  multiple- star  formation. 
If  chaotic  Y-body  dynamics  is  the  dominant  mechanism,  sub¬ 
systems  in  the  secondaries  should  be  much  less  common  than 
in  the  primaries  because,  in  triple  encounters,  the  lightest  star 
is  usually  ejected  while  the  two  most  massive  stars  pair  in  a  bi¬ 
nary.  On  the  other  hand,  the  hydrodynamic  simulations  of  Bate 
(2012)  predict  a  large  fraction  of  multiple  systems,  including 
2+2  quadruples  (i.e.,  two  close  pairs  in  a  wide  outer  binary). 

Recent  advances  in  observational  techniques  have  substan¬ 
tially  improved  our  knowledge  of  the  distribution  of  binary  pe¬ 
riods  and  mass  ratios,  as  reviewed  by  Duchene  &  Kraus  (2013). 
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Robo-AO  Observation  Parameters 
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Figure  1.  Cumulative  histograms  of  the  Ic  magnitudes  of  components  observed 
with  Robo-AO. 


Lucky  imaging  has  been  used  to  survey  large  samples  (Ginski 
et  al.  2012;  Janson  et  al.  2012;  Jodar  et  al.  2013).  Yet,  the  em¬ 
phasis  is  almost  always  placed  on  binaries,  leaving  higher-order 
hierarchies  aside  or  addressing  them  casually;  a  rare  exception 
is  the  work  of  Law  et  al.  (2010)  who  observed  a  higher  than 
expected  high-order  multiplicity  of  M  dwarf  stars. 

This  survey  aims  at  complementing  the  statistics  of  hierar¬ 
chical  systems  in  the  FG-67  sample  by  imaging  with  Robo-AO 
(Baranec  et  al.  2014).  Our  main  goal  is  to  constrain  the  mul¬ 
tiplicity  of  secondary  components,  most  of  which  so  far  had 
no  high-resolution  imaging  data.  The  second  goal  is  to  search 
for  additional  companions  around  primary  targets  that  are  them¬ 
selves  known  close  binaries.  Here  we  explore  the  poorly  covered 
part  of  the  parameter  space  at  separations  on  the  order  of  an  arc- 
second,  searching  for  low-mass  tertiary  components.  In  short, 
we  focus  on  binaries  (both  wide  and  close)  and  quantify  the 
frequency  of  additional  companions. 

We  begin  with  a  brief  characterization  of  the  surveyed  sample 
in  Section  2.  Observations  with  Robo-AO  and  data  processing 
are  covered  in  Section  3,  with  data  tables  described  in  Section  4. 
In  Section  5  the  results  of  this  survey  are  presented,  and  Section  6 
is  the  summary. 

2.  THE  TARGET  SAMPLE 

The  targets  for  this  survey  were  selected  from  the  FG-67 
database  (Tokovinin  2014).  We  observed  known  and  suspected 
binary  stars  and  tried  to  constrain  the  frequency  of  additional 
components  in  those  systems;  high  priority  was  placed  on  distant 
secondary  components  in  wide  binaries.  The  selection  criteria 
were  separations  of  >10''  and  a  declination  north  of  —15°. 
Some  brighter  secondaries  are  in  the  Hipparcos  catalog,  and 
therefore  “screened”  for  companions  by  that  space  mission  (but 
not  as  deeply  as  with  Robo-AO).  The  majority,  however,  are 
rather  faint.  Many  of  those  secondaries  were  discovered  only 
recently  by  dedicated  programs  (Tokovinin  2011;  Tokovinin  & 
Lepine  2012). 

The  second-priority  list  contains  mainly  targets  that  are  bi¬ 
naries  with  periods  P  <  100  yr  (separation  ~0'.'5  at  50  pc 
distance).  The  rationale  for  restricting  the  periods  is  related 
to  dynamical  stability:  the  semi-major  axis  of  the  outer  orbit 
should  be  at  least  three  times  larger  than  that  of  the  inner  or¬ 
bit,  and  should  match  well  the  discovery  space  of  Robo-AO. 
Harrington  (1992)  quantified  a  value  of  three  for  the  ratio  of 


Observation  date  range 

2012  Aug  31  to  2013  August  21 

Targets  observed 

695 

Total  observing  time 

29  hr 

Observation  filters 

SDSS  i',  LP600 

Exposure  time 

60s 

Typical  FWHM  resolution 

0'.'12to0'.'15 

Typical  contrast  ratio  at  O'.'  5 

A  ~  5  mag. 

periastron  of  the  outer  binary  to  apastron  of  the  inner  binary 
as  the  critical  factor  for  long-term  stability  (assuming  equal 
masses).  The  program  was  later  complemented  with  observa¬ 
tions  of  additional  resolved  visual  binaries  with  separations 
from  0'.'2  to  2",  and  P  >  100  yr.  Here  the  chances  of  dis¬ 
covering  wider  tertiary  components  are  small.  The  chances  of 
resolving  a  secondary  into  a  hereto  unknown  close  pair  are  also 
small.  However,  observing  those  stars  (mostly  resolved)  was 
useful  for  improving  their  orbits  and  for  confirming  some  un¬ 
certain  visual  pairs  that  were  resolved  only  once  in  the  past  and 
never  confirmed  since.  The  program  included  repeated  obser¬ 
vations  of  several  known  binaries  as  calibrators  for  data  quality 
control. 

We  observed  most,  but  not  all,  components  from  the  original 
lists:  239  secondaries  (212  of  those  with  good  quality),  354  close 
pairs,  and  102  resolved  binaries.  Some  secondaries  turned  out  to 
be  too  faint.  Figure  1  presents  the  distribution  of  Iq  (Cousins  7) 
magnitudes  of  our  primary  targets  (van  Leeuwen  2007)  and  of 
their  secondary  components.  The  latter  are  derived  from  their 
J  magnitudes  in  the  Two  Micron  All  Sky  Survey  (2MASS) 
catalog  (Cutri  et  al.  2003),  assuming  that  the  stars  are  on  the  main 
sequence.  All  secondaries  with  useful  data  have  Ic  <  13.5  mag. 
The  median  masses  of  primary  and  secondary  components 
estimated  from  their  absolute  magnitudes  are  1 .22  and  0.7 1  A4o , 
respectively;  80%  of  secondary  masses  are  between  0.31  and 
1.04  Mo. 

3.  OBSERVATIONS  AND  DATA  REDUCTION 

The  instrument  used  for  this  survey  was  Robo-AO,  the  first 
robotic  laser  guide  star  AO  system  (Baranec  et  al.  2014), 
which  is  designed  to  operate  automatically  on  1-3  m  class  tele¬ 
scopes  in  order  to  undertake  high  efficiency  observing  programs 
(e.g.,  large  surveys).  The  prototype  Robo-AO  instrument  is  cur¬ 
rently  deployed  on  the  Palomar  Observatory  60  inch  (1.5  m) 
telescope  (P60). 

As  a  fully  automated  system,  Robo-AO  is  unique  in  its  ability 
to  observe  targets  at  both  high  resolution  and  high  cadence. 
Robo-AO  uses  Rayleigh  scattering  from  a  UV  laser  focused 
at  10  km  from  the  telescope  as  the  wavefront  reference,  and 
generates  images  at  the  diffraction  limit  of  the  P60  (O'.' 12-0'.' 15 
FWHM)  with  Strehl  ratios  of  10%-25%  in  the  visible  filters 
used  for  this  survey.  The  AO  system  corrects  the  high  order 
wavefront  aberrations  with  automated  software  that  operates  at 
a  rate  of  1.2  kHz,  sharpening  the  instantaneous  point-spread 
function  (PSF)  across  the  science  camera  field  of  view  (FOV). 
A  bright  star  within  the  FOV  is  still  required  to  correct  the  tip- 
tilt  motion;  this  is  achieved  in  the  automated  data  processing 
software  (Section  3.2).  Typically,  Robo-AO  requires  a  star 
brighter  than  m  ~  16  mag  with  a  broadband  filter  to  apply 
the  tip-tilt  correction  successfully.  Instrument  parameters  for 
this  survey  are  listed  in  Table  1 . 
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3.1.  Observations 

Baranec  et  al.  (2013)  and  Baranec  et  al.  (2014)  describe 
the  operation  and  the  science  instrumentation  of  the  Robo-AO 
system.  The  output  of  the  science  camera  raw  data  are  image 
cubes  composed  of  1024  x  1024  pixel  image  frames;  a  total 
of  516  frames  were  gathered  during  each  60  s  exposure,  which 
were  then  combined  into  a  single  image  for  further  analysis 
by  the  automated  data  processing  software.  Almost  all  targets 
were  observed  with  the  Sloan  Digital  Sky  Survey  (SDSS)  i' 
filter  (York  et  al.  2000)  with  a  central  wavelength  of  754  nm 
and  a  FWHM  bandpass  of  119  nm.  Some  fainter  targets  were 
observed  with  a  long  pass  600  nm  filter  (LP600)  that  transmits 
all  wavelengths  longer  than  600  nm  to  the  quantum  efficiency 
cutoff  of  the  detector. 

The  pointing  error  of  the  P60  telescope  with  Robo-AO 
mounted  is  on  the  order  of  10'',  and  can  vary  between  observing 
runs  depending  on  telescope  balance  and  adjustments.  With  a 
44"  FOV,  and  no  clear  way  to  select  one  target  over  another 
in  the  stellar  field,  Robo-AO  does  not  attempt  to  recenter  the 
target  automatically  in  fear  of  selecting  the  wrong  star  and 
moving  the  science  target  out  of  the  field.  This  effect  causes 
a  large  variation  in  the  placement  of  the  target  star  in  the 
image  frame;  a  planned  upgrade  to  the  Robo-AO  software  and 
continuing  improvements  to  the  P60  will  minimize  this  effect. 
All  observations  were  checked  to  confirm  that  the  science  target 
was  captured  in  the  CCD  frame.  The  consequence  for  this  survey 
is  that  target  placement  is  across  the  entire  FOV,  and  targets 
sometimes  ended  up  uncomfortably  close  to  the  edge  of  the 
frame,  limiting  the  observable  area  around  the  target.  When 
necessary,  stars  affected  were  reobserved  to  capture  a  larger 
area  and  allow  a  better  examination  for  companions. 

Robo-AO  can  observe  targets  at  a  rate  of  "^20  per  hour, 
with  an  intelligent  queue  system  selecting  the  best  object  to 
observe  at  a  given  time  (Riddle  et  al.  2014).  The  Robo-AO 
queue  system  interleaves  several  different  science  programs 
through  a  single  night.  Observations  of  this  survey  were  spread 
over  a  year  (2012  July  to  2013  August)  to  gather  targets  at 
all  hours  of  right  ascension.  The  entire  695  targets  of  this 
survey  used  29  total  hours  of  observing  time  (this  time  also 
included  calibration  binary  observations  and  reobservations). 
Robo-AO  is  ideally  suited  to  a  large  survey  such  as  this  one, 
and  is  currently  the  only  AO  system  that  can  observe  this 
many  targets  in  such  a  short  time  (Terziev  et  al.  2013;  Law 
et  al.  2014a). 

3.2.  Automated  Data  Processing 

The  data  were  processed  by  automated  reduction  software 
developed  for  Robo-AO  image  alignment  (see  Law  et  al.  2014b, 
2014a  for  the  details).  Each  individual  frame  is  corrected  for 
bias,  dark  current  and  fiat  field  using  standard  calibration  data, 
and  the  frames  are  then  run  through  the  software  to  co-add 
them  and  create  a  final  high  resolution  image  oversampled  to 
2048  X  2048. 

Stars  in  each  output  image  were  identified  visually  and  their 
approximate  centers  were  marked.  Relative  astrometry  and  pho¬ 
tometry  of  wide  pairs  where  the  images  do  not  overlap  was  done 
by  fitting  the  scaled  and  shifted  image  of  the  brightest  compo¬ 
nent  to  the  secondary  within  10  pixel  radius  from  the  image 
center.  For  closer  (overlapping)  binaries  we  used  blind  decon¬ 
volution,  as  described  by  Tokovinin  et  al.  (2013).  However, 
most  of  these  “blind”  measures  are  superseded  by  the  results  of 
speckle  processing. 


3.3.  Binary -star  Measurement  by  Speckle  Processing 

Speckle  processing  is  complementary  to  the  blind  deconvo¬ 
lution  of  images,  as  it  delivers  diffraction-limited  resolution 
even  for  low-Strehl  data,  when  the  instantaneous  PSF  has  mul¬ 
tiple  spikes  (speckles).  The  algorithm  takes  care  of  the  multi¬ 
speckle  PSF  and  works  well  even  at  low  flux,  when  the  selection 
of  the  brightest  pixel  for  recentering  lucky  images  is  compro¬ 
mised  by  the  photon  noise.  Moreover,  the  power  spectrum  (or 
auto-correlation  function,  ACF)  is  proportional  to  the  square  of 
the  signal  and,  therefore,  it  automatically  assigns  high  weight 
to  sharp  images.  For  widely  spaced  binaries,  speckle  process¬ 
ing  is  done  for  each  component  individually  to  detect  close 
subsystems. 

Subframes  of  256  x  256  original  (unbinned)  pixels  cen¬ 
tered  on  each  component  were  selected  from  the  data  cubes. 
The  background  was  estimated  as  median  signal  in  each  col¬ 
umn  of  the  full  image  and  subtracted,  removing  the  systematic 
bias  pattern  along  the  CCD  lines.  No  fiat-field  correction  was 
applied.  The  average  power  spectrum  was  calculated  on  the 
256  frame  cube  for  each  selected  component.  In  parallel,  an  im¬ 
age  was  produced  from  these  data  by  recentering  on  the  selected 
component  and  weighting  individual  frames  with  the  maximum 
intensity  in  each  frame.  The  weight  is  therefore  proportional  to 
the  image  sharpness.  These  auxiliary  images  are  not  used  for 
astrometry  or  photometry,  but  are  helpful  for  verifying  compan¬ 
ion  detection  and  for  resolving  the  180°  ambiguity  of  position 
angle  inherent  to  the  classical  speckle  processing. 

The  algorithm  of  speckle  data  processing  and  binary-star 
measurement  used  here  is  adapted  from  the  work  of  Tokovinin 
et  al.  (2010).  Figure  2  illustrates  the  case  of  a  close  0'.'12  pair 
near  the  diffraction  limit.  Extracting  astrometry  and  photometry 
from  the  image  appears  problematic,  whereas  the  “fringes”  in 
the  power  spectrum  are  very  clear  and  constrain  the  position  and 
magnitude  difference  Am. 

For  binaries  with  nearly  equal  components,  the  pipeline  image 
contains  an  “antipode”  because  the  brightest  pixel  on  which  the 
frames  are  recentered  may  belong  to  either  binary  component. 
In  this  case  the  relative  photometry  derived  from  the  image 
is  wrong  (although  correctable),  but  the  speckle  processing 
delivers  the  correct  Am .  It  also  increases  the  resolution  on  distant 
secondary  components  where  the  full-frame  images  suffer  from 
the  tilt  anisoplanatism.  On  the  other  hand,  if  the  secondary 
component  is  very  faint,  the  pipeline  technique  of  recentering 
on  the  bright  primary  works  much  better  than  speckle  processing 
of  the  secondary  alone.  We  quantify  the  strength  of  the  speckle 
signal  by  its  ratio  to  the  level  of  photon  noise  at  a  spatial 
frequency  two  times  less  than  the  cutoff  frequency,  and  consider 
all  measures  resulting  from  the  “weak”  data  to  be  uncertain. 

Speckle  processing  and  imaging  are  truly  complementary. 
In  the  final  step,  we  combine  binary-star  measures  from  those 
two  processing  techniques.  For  binaries  with  p  <  2"  the  speckle 
measures  are  preferred  (with  a  few  exceptions  where  the  speckle 
processing  failed),  for  wider  pairs  measures  from  the  images 
are  retained.  Comparison  of  repeated  measures  of  the  same 
pairs  by  speckle  and  blind  deconvolution  shows  their  excellent 
agreement  and  gives  an  idea  of  the  precision.  Table  2  lists 
average  values  and  rms  scatter  for  two  sample  binaries  measured 
by  both  methods  several  times. 

3.4.  Detection  Limits 

Limits  for  detecting  companions  in  the  pipeline  images  are 
evaluated  in  the  standard  way.  Fluctuation  of  the  signal  in 
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Table  2 

Comparison  between  Speckle  and  Blind  Deconvolution 


HIPl 

Discoverer 

Xobs 

Speckle  (mean,  rms) 

Blind  (mean,  rms) 

6 

P 

Am 

6 

P 

Am 

20472 

STF  535 

9 

271.15 

1.130 

0.98 

271.18 

1.129 

1.04 

0.32 

0.002 

0.10 

0.30 

0.002 

0.06 

23100 

COU  1525 

8 

293.55 

1.833 

3.36 

293.56 

1.832 

3.32 

0.21 

0.007 

0.05 

0.24 

0.007 

0.12 

annular  zones  surrounding  the  star  are  computed  and  it  is 
assumed  that  the  companions  brighter  than  5a  are  detectable. 
The  same  method  is  applied  to  the  ACFs  in  the  speckle 
processing.  The  two  estimates  of  the  detection  limits  obtained 
in  this  way  are  very  similar,  with  the  speckle  limits  being 
normally  a  little  deeper  (Figure  3).  Depending  on  variable 
AO  correction  and  target  brightness,  the  individual  detection 
limits  vary  substantially.  For  each  target,  we  list  the  best 
(deepest)  limits  at  three  characteristic  separations  of  O'.' 15, 
O'.'S  and  2'.'1,  selecting  the  best  of  two  alternative  processing 
techniques.  When  the  target  was  observed  several  times,  the 
deepest  detection  limits  are  reported  in  the  table  of  final  results. 

All  of  the  survey  primary  targets  are  bright,  but  some 
secondary  components  present  a  problem  as  they  are  too  faint 
to  provide  a  useful  signal  in  the  power  spectrum.  Their  pipeline 
images  can  show  a  bright  1  pixel  spike  at  the  center  surrounded 
by  a  fuzzy  halo;  this  is  created  by  recentering  on  the  photon 
events,  rather  than  on  the  brightest  speckle,  so  the  resolution 
is  lost.  Detection  limits  computed  formally  from  such  images 
are  over-estimated.  However,  we  prefer  to  keep  such  weak  data 
in  the  final  table,  despite  this  obvious  caveat,  because  they  still 
contain  some  information  on  the  duplicity. 


3.5.  Calibration  and  Distortion  Correction 

The  reimaging  system  of  the  Robo-AO  instrument  contains  a 
double  optical  relay  with  off-axis  parabolic  mirrors.  Such  relays 
are  known  to  introduce  quadratic  image  distortion.  In  the  case 
of  Robo-AO,  the  distortion,  as  reported  by  the  optical  design, 
is  quite  strong.  It  is  directed  along  the  visible  CCD  columns 
and  displaces  all  sources  down  by  as  much  as  26  pixels  in  the 
corners. 
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Figure  3.  Median  detection  limit  Am(p)  determined  on  the  images  (full  line) 
and  on  the  speckle  ACFs  (dashed  line).  Measured  pairs  are  plotted  as  crosses. 

Independently  of  the  optical  design,  the  distortion  in 
Robo-AO  was  mapped  on  the  sky  using  the  image  of  the  globular 
cluster  M  15  (S.  R.  Hildebrandt,  private  communication,  2012). 
The  pixel  scale  in  the  X-direction  was  found  to  be  43.74  mas, 
while  the  scale  in  Y  was  slightly  different.  The  orientation  is 
such  that  the  -fF  axis  of  the  detector  points  at  a  position  an¬ 
gle  23?  9  from  North  and  the  +X  axis  points  East  (the  image 
has  mirror  orientation).  The  quadratic  distortion  agrees  with  the 
optical-design  data  and  is  directed  along  Y. 

If  we  express  the  actual  A,  F  star  coordinates  in  the  image 
in  the  “small”  over-sampled  pixels  of  the  automatic  pipeline 
reduced  image,  the  undistorted  coordinate  F'  is  determined  as 

Y'  =  Y  -  0.0299Ay  +  Cxi^Xf  +  CyiAYf,  (1) 

where  AX  and  AF  are  counted  from  the  pipeline  image  center 
at  (1024,1024)  and  the  coefficients  are  Cy  =  —2.54  10“^, 
Cx  =  —2.55  10“^.  The  linear  term  in  F  corresponds  to  the 
modified  vertical  scale.  It  is  caused  by  the  mismatch  between  the 
optical  and  CCD  centers  combined  with  the  quadratic  distortion. 

This  distortion  substantially  affects  the  measurements  of  even 
close  binaries  because  the  CCD  lines  projected  on  the  sky 
are  parabolic.  The  distortion-induced  tilt  of  the  lines  reaches 
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Table  3 

List  of  Observed  Components  and  Detection  Limits 


HIPl 

Comp. 

RA  (2000) 

m 

Dec  (2000) 

(0  / 

y 

(mag) 

Sep. 

(") 

Epoch 

+2000 

a'i5 

(mag) 

a'8 

(mag) 

2'.T 

(mag) 

Pmax 

(") 

179 

A 

00  02  16.68 

-13  24  26.8 

6.90 

0.0 

13.5730 

2.7 

5.2 

6.3 

13.9 

223 

A 

00  02  47.17 

+02  07  48.5 

7.17 

0.0 

12.7049 

2.1 

4.9 

5.7 

17.5 

394 

A 

00  04  53.76 

+34  39  35.3 

6.09 

0.0 

12.7020 

1.7 

5.1 

6.7 

15.1 

493 

D 

00  05  54.75 

+18  14  06.0 

7.47 

0.0 

12.6747 

1.2 

4.5 

6.9 

11.8 

601 

A 

00  07  18.15 

+20  57  54.3 

9.14 

0.0 

12.7050 

1.8 

4.2 

4.7 

17.9 

754 

B 

00  09  13.92 

+25  16  38.9 

11.46 

29.6 

12.6747 

1.6 

4.7 

6.8 

12.2 

(This  table  is  available  in  its  entirety  in  machine-readable  form.) 


5? 6  in  the  corners  of  the  CCD.  We  correct  the  geometric 
distortion  using  the  pixel  coordinates  of  each  source  in  the 
images.  Equation  1  is  applied  to  each  component  of  a  binary. 
The  parameters  p,0  are  recalculated  from  the  difference  of 
undistorted  pixel  coordinates  X,  Y'.  In  some  instances  the  angles 
of  close  binaries  are  changed  by  as  much  as  2°  after  distortion 
correction.  The  scatter  of  relative  positions  of  binaries  measured 
several  times  is  reduced  after  distortion  correction,  while  the 
relative  positions  of  wide  pairs  became  closer  to  their  positions 
derived  from  2MASS. 

3.6.  Caveats 

As  with  any  observing  program,  there  were  issues  with  some 
of  the  data  collected  that  required  some  extra  effort  to  make  it 
useful.  Each  pipeline  image  was  examined  manually  to  find  and 
remove  erroneous  measurements  or  detections;  a  graphical  IDE 
tool  was  created  to  do  this,  and  also  used  to  mark  the  components 
and  fit  wide  binaries.  Owing  to  the  P60  pointing  errors,  some 
images  were  empty,  or  the  targets  were  found  close  to  the  edge 
of  the  EOV  (see  Section  3.1). 

The  quality  of  AO  compensation,  Strehl  ratio,  and  the 
width  of  the  PSE  depended  on  the  seeing  conditions  and 
varied  substantially.  As  a  consequence,  the  depth  of  companion 
detection  was  also  variable.  In  some  instances  we  detected  new 
companions  in  the  good-quality  images,  but  missed  them  in  the 
poor  ones.  The  PSE  had  a  persistent  structure  (“static  speckle”) 
that  could  be  mistaken  for  a  companion.  This  structure  also 
depended  on  the  AO  compensation  quality  and,  possibly,  on 
other  factors.  Eortunately,  there  were  always  other  stars  observed 
in  the  same  conditions  so  the  PSE  could  be  compared  in  order 
to  verify  suspicious  detections. 

In  a  few  cases,  the  pipeline  algorithm  occasionally  centered 
a  frame  on  a  bright  spike  caused  by  cosmic  rays,  while  the 
remaining  frames  were  centered  correctly.  This  produced  a  false 
satellite  that  could  be  mistaken  for  a  new  companion.  Also, 
several  faint  satellites  caused  by  internal  reflections  in  the  optics 
are  always  visible  around  bright  (V  <  3  mag)  targets.  These 
reflections  effectively  reduce  the  dynamic  range  of  companion 
detection.  We  were  able  to  reduce  some  images  where  the 
primary  component  was  saturated.  Such  data  are  still  useful 
for  detecting  wide  faint  companions.  Relative  astrometry  and 
photometry  of  such  binaries  was  done  by  fitting  PSE  within 
an  annulus,  i.e.,  excluding  the  central  saturated  PSE  core  from 
the  fit. 

The  accuracy  of  the  measured  positions  of  wide  binaries  is 
less  than  for  the  close  ones.  This  is  likely  caused  by  the  residual 
tilt  anisoplanatism:  the  differential  tilt  between  the  components 
caused  by  high-altitude  turbulence  is  not  completely  averaged 
during  the  short  60  s  exposure. 


In  many  cases,  we  were  able  to  reobserve  targets  after  the 
data  analysis  showed  that  the  data  were  not  sufficient;  the  ability 
of  Robo-AO  to  observe  quickly  and  efficiently  allowed  many 
objects  to  be  added  to  the  final  analysis  that  would  have  not  been 
available  otherwise. 

4.  OBSERVING  DATA  TABLES 

The  list  of  all  components  successfully  observed  in  this  pro¬ 
gram  is  given  in  Table  3.  The  list  contains  695  components 
belonging  to  595  systems.  Its  columns  contain  the  HIPl  designa¬ 
tion  (which  is  the  Hipparcos  number  of  the  primary  component 
in  each  system,  and  is  used  as  the  main  identifier  through  this 
paper),  component  designation,  equatorial  coordinates  of  the 
observed  component,  its  V  magnitude,  and  its  separation  from 
the  primary  (these  data  are  taken  from  the  main  database  of  the 
EG-67  sample).  Secondary  components  have  non-zero  separa¬ 
tions.  The  following  columns  list  the  date  of  observation  and 
detection  limits  at  separations  of  O'.'! 5,  WS,  and  2('l.  The  last 
column  gives  the  distance  from  the  star  to  the  nearest  frame 
border,  indicating  the  maximum  separation  of  its  detectable 
satellites. 

Measurements  of  resolved  doubles  are  given  in  Table  4.  Each 
pair  is  identified  by  its  WDS  code  (Mason  et  al.  2001b),  dis¬ 
coverer  designation,  and  component  designation.  In  column  3, 
the  HIPl  number  of  the  primary  star  in  each  system  is  given. 
Eollowing  that  are  the  date  and  number  of  observations  (multi¬ 
ple  observations  made  in  the  same  filter  within  0.2  yr  are  aver¬ 
aged),  position  angle  (degrees),  position  angle  error  (degrees), 
separation  (arcseconds),  separation  error(milliarcseconds),  and 
magnitude  difference.  The  errors  of  position  angle  and  separa¬ 
tion  are  computed  as  errors  of  the  mean  in  the  case  of  averaged 
data,  otherwise  estimated  from  the  speckle  processing,  or  listed 
as  zero  for  PSP  fitting  of  wide  pairs.  Uncertain  measures  (weak 
speckle  signal  or  blind  deconvolution)  are  marked  by  colons 
after  Am.  We  also  mark  the  bandpass  when  it  differs  from  V 
(the  LP600  filter  is  denoted  by  w).  The  last  three  columns  give 
deviations  in  position  angle  and  separation  from  the  orbits,  when 
available. 

5.  RESULTS 

5.1.  Newly  Resolved  Systems 

The  list  of  newly  resolved  binary  companions  (including 
unrelated  background  stars)  is  extracted  from  the  main  table  and 
presented  separately  in  Table  5;  they  are  given  the  “discoverer 
code”  RAO  in  the  WDS;  a  space  is  added  between  discoverer 
codes  and  component  designations  to  make  it  clear  they  are 
not  the  same  (i.e.,  component  designations  such  as  Aa  and 
discoverer  codes  such  as  RAO  20  are  parallel  and  must  not 
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Figure  4.  Images  of  some  new  binaries  discovered  with  Robo-AO.  Each  fragment  of  the  pipeline  image  has  a  size  of  100  pixels  (272),  the  intensity  scale  is  adjusted 
to  highlight  the  companions,  indicated  by  arrows.  The  images  are  displayed  without  rotation  and  labeled  by  HIP  numbers. 
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Figure  5.  Power  spectra  of  some  newly  resolved  close  binaries.  Primaries  are  identified  by  their  HIPl  number. 


Table  4 

Measurements  of  Binary  Stars  with  Robo-AO 


WDS 

(2000) 

Discoverer 

Designation 

HIPl 

Epoch 

+2000 

N 

e 

(deg) 

ere 

(deg) 

p 

(") 

ap 

(mas) 

Am 

(mag) 

[0-C]e 

(deg) 

[0-C|p 

(") 

Ref. 

00028+0208 

BU281  AB 

223 

12.7049 

1 

160.8 

0.0 

1.567 

0.5 

1.66 

00057+4549 

STT547AB 

473 

12.7723 

3 

186.9 

0.2 

6.026 

42.8 

0.09 

-0.3 

0.001 

1 

13.0478 

2 

186.4 

0.1 

5.758 

68.7 

0.16 

-0.9 

-0.266 

1 

00057+4549 

POP  217  AP 

473 

13.0478 

2 

325.5 

0.3 

13.140 

49.4 

5.03 

00073+2058 

HDS12 

601 

12.7050 

1 

6.7 

0.1 

1.864 

2.0 

2.38: 

00251+4803 

HDS56AB 

1987 

12.7020 

1 

279.3 

0.1 

0.324 

0.5 

2.19 

-0.8 

-0.016 

2 

00261-1123 

YR4 

2066 

13.5730 

1 

176.8 

0.0 

0.376 

0.0 

3.12: 

00271-0753 

A431  AB 

2143 

13.5730 

1 

157.6 

0.1 

0.166 

0.3 

0.00 

-5.0 

-0.019 

3 

References.  (1)  Kiyaeva  et  al.  (2001);  (2)  Cvetkovic  (2011);  (3)  Scardia  (1981). 
(This  table  is  available  in  its  entirety  in  machine-readable  form.) 


be  confused).  Some  binaries  were  measured  several  times, 
confirming  the  measured  resolution.  Images  of  new  close  pairs 
are  presented  in  Figure  4.  The  closest  pairs  are  detected  in 
the  speckle  processing  by  examining  the  power  spectra;  these 


cases  are  illustrated  by  Figure  5.  The  new  wide  physical 
pairs  are  confirmed  by  comparing  the  measured  positions  of 
the  components  with  their  position  in  the  2MASS  catalog, 
and  by  considering  additional  information  such  as  magnitude. 
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Table  5 

New  Pairs  Resolved  with  Robo-AO 


WDS 

(2000) 

HIPl 

Discoverer 

Designation 

N 

0 

(deg) 

p 

(") 

Am 

(mag) 

Status,  remarks 

00293-0555 

2292 

RAO  1 BC 

5 

166.7 

0.53 

3.3 

P  (fixed  for  1  yr),  B  =  HIP  2350 

00346+6235 

2717 

RA02BC 

1 

612 

3.55 

4.3 

P  (2MASS  image) 

00346+6235 

2717 

RA02BD 

1 

302.5 

3.180 

6.85 

O  (crowded  field,  non-hierarchical) 

00409+3107 

3214 

RA0  3Ba,Bb 

3 

275.6 

1.629 

0.0 

P  (bright,  2MASS  image,  fixed  for  1  yr) 

00452+0015 

3540 

RA04BC 

1 

223.1 

1.521 

3.5 

P  (confirmed  with  P200,  low  crowding) 

00487+1841 

3795 

RA05Ca,Cb 

1 

322.4 

0.674 

1.3 

P  (bright  &  close,  low  crowding) 

00515+5630 

4016 

RA0  38BC 

1 

270.4 

16.999 

3.97 

O  (crowded  field,  2MASS  position) 

00515+5630 

4016 

RAO  38  CD 

1 

82.5 

1.720 

1.28 

O?  CD  is  unrelated  to  B 

01027+0908 

4878 

RA0  39AB 

1 

36.2 

2.927 

6.57 

O?  (different  position  in  1  yr) 

01027+0908 

4878 

RAO  39  AC 

1 

147.0 

32.155 

6.33 

O  (2MASS  position,  small  PM) 

01075+4116 

5276 

RAO40AB 

2 

335.4 

6.167 

6.4 

P  (2MASS  position,  conf.  with  P200) 

01080+3915 

5313 

RA041  AB 

1 

244.7 

26.24 

4.4 

P  (2MASS  position,  color) 

01253+0128 

6653 

RA06AB 

3 

204.5 

1.827 

4.2 

P  (low  crowding,  red  color) 

01409+6410 

7845 

RAO  32  AC 

5 

175.3 

24.84 

6.5 

O  (fast  relative  motion),  not  in  FG-67 

02308+5533 

11696 

RAO  42  BE 

1 

207.0 

15.677 

5.11 

O  (2MASS  position,  high  crowding) 

02355+5948 

12062 

RA0  43AB 

1 

207.8 

14.472 

6.15 

O  (2MASS  position,  high  crowding) 

02356+6106 

12067 

RA07AB 

1 

169.0 

5.71 

5.1 

P  (2MASS  position) 

02370+2439 

12189 

RA08BC 

2 

285.2 

0.536 

4.3 

P  (fixed  position,  conf.  with  P200).  B  =  HIP  12184 

02430+5812 

12685 

RA044AB 

1 

211.1 

21.670 

4.08 

O?  (color,  high  crowding) 

02462+0536 

12925 

RA09AB 

3 

252.8 

1.91 

4.4 

P  (fixed  position,  low  crowding) 

02517+4559 

13336 

RAO  45  AC 

1 

81.8 

15.195 

6.30 

O  (2MASS  position) 

03177+3838 

15329 

RAO10Ca,Cb 

3 

223.4 

1.560 

4.0 

P  (small  separation,  red  color) 

03305+2006 

16329 

RA011Ba,Bb 

4 

75.9 

0.159 

0.9 

P  (spectroscopic  binary),  B  =  HD  21663 

03390+4232 

17022 

RA047AB 

1 

142.9 

1.741 

6.2 

P  (confirmed  with  P200) 

03413+4554 

17217 

RAO  48  AC 

1 

269.0 

4.83 

6.4 

P  (2MASS  image) 

03426+1718 

17336 

RA078AB 

1 

269.7 

10.388 

6.44 

O  (2MASS  position) 

03537+5316 

18218 

RAO  49  AC 

1 

241.6 

25.912 

6.65 

O  (2MASS  position,  high  crowding) 

04007+2023 

18719 

RAO50AB 

1 

308.4 

5.905 

6.35 

O  (not  found  in  2MASS  images) 

04092+4010 

19389 

RAO  79  CD 

2 

294.3 

9.584 

3.04 

O  (2MASS  position) 

04313+2008 

21099 

RA051  AB 

1 

310.3 

18.800 

6.32 

O  (2MASS  position) 

04313+2008 

21099 

RA051BC 

1 

177.4 

0.917 

0.26 

O?  (unrelated  to  A) 

04363+5502 

21443 

RA0  35AB 

1 

14.4 

5.737 

4.65 

P?  (2MASS  position  &  color,  confirmed  with  P200) 

05017+2050 

23396 

RAO  52  AC 

1 

242.8 

38.414 

3.85 

O  (2MASS  position) 

05096+2947 

24016 

RA053BC 

1 

106.2 

26.153 

6.58 

O  (2MASS  position),  B  =  HIP  24005 

05247+6323 

25300 

RAO  36  AC 

1 

227.1 

6.84 

5.5 

P  (2MASS  image,  conf.  with  P200) 

05376+0607 

26444 

RA054BC 

1 

34.9 

4.526 

5.79 

O  (not  in  2MASS  image) 

05376+0607 

26444 

RA054BD 

1 

157.2 

4.398 

6.12 

O  (not  in  2MASS  image) 

05444+4024 

27067 

RAO  55  AC 

1 

151.9 

11.103 

5.55 

O  (2MASS  position,  high  crowding) 

05465+7437 

27246 

RAO  37  AC 

1 

68.5 

10.85 

3.5 

P  (2MASS  position,  color) 

06335+4822 

31267 

RAO  80  AC 

1 

221.9 

4.97 

5.3 

P  (2MASS  image,  conf.  with  P200) 

06562+4032 

33355 

RA056AB 

1 

156.6 

5.510 

6.35 

P?  (2MASS  image,  conf.  with  P200) 

07043-0303 

34110 

RAO  57  AC 

1 

10.6 

15.713 

6.64 

O  (2MASS  position) 

07172+3306 

35265 

RA058BC 

1 

126.6 

17.637 

5.45 

O  (2MASS  position,  high  crowding) 

07475+7012 

38018 

RAO  59  AC 

1 

118.8 

15.476 

5.87 

O  (2MASS  position) 

08138+6306 

40298 

RAO60BC 

2 

37.6 

13.074 

3.48 

O  (2MASS  position) 

08159+0227 

40479 

RA061  AB 

1 

170.6 

31.245 

6.42 

O  (2MASS  position) 

08211+6527 

40918 

RA013BC 

1 

245.6 

2.79 

3.2 

P  (2MASS  image),  B  =  HIP  40882 

08258+1703 

41319 

RA062BC 

2 

310.0 

16.798 

4.19 

O  (2MASS  position) 

08508+3504 

43426 

RAO  63  CD 

1 

357.6 

18.943 

3.20 

O  (2MASS  position) 

09505+0421 

48273 

RAO90Ba,Bb 

1 

83.8 

0.160 

0.0 

P  (close  pair) 

10079+6650 

49638 

RAO  15  AB 

1 

71.4 

0.510 

3.0: 

P  (close  pair,  low  crowding) 

14094+1015 

69160 

RA016AB 

1 

261.2 

0.119 

0.7 

P  (SB 2  resolved) 

14113+5424 

69322 

RA017BC 

2 

145.9 

0.577 

1.7 

P  (close  pair,  low  crowding) 

14418+2432 

71843 

RA064BC 

2 

61.9 

17.271 

2.02 

O  (blue  color) 

15277+4253 

75676 

RA018BC 

1 

100.2 

0.408 

1.5 

P  (close  pair,  small  Am) 

16147+3352 

79607 

RAO  65  EF 

1 

238.2 

32.267 

3.74 

O  (not  found  in  2MASS),  E  =  HIP  79551 

16150+6040 

79629 

RA0  81Ba,Bb 

1 

42.1 

0.139 

0.0 

P  (close  pair,  small  Am) 

16400-0605 

81608 

RA0  82BC 

1 

148.9 

27.160 

3.89 

O  (2MASS  position) 

17229-0223 

85042 

RA019BC 

1 

86.3 

0.752 

0.0 

P  (close  pair,  small  Am) 

17422+3804 

86642 

RAO20AB 

1 

302.2 

2.20 

5.0 

P  (confirmed  with  P200) 

18123+0154 

89207 

RAO  66  AC 

2 

71.9 

17.237 

6.81 

O?  (2MASS  position,  crowded  field) 

18352+4135 

91120 

RA0  83Ba,Bb 

1 

336.1 

0.136 

0.8 

P  (close  pair,  small  Am) 

19145+3434 

94540 

RA0  84AE 

1 

72.2 

2.900 

4.50 

O?  (not  confirmed  in  2MASS  image) 

19158+3823 

94666 

RAO  85  AB 

1 

331.1 

3.57 

5.7 

P  (2MASS  image) 

19188+1629 

94905 

RA067AB 

1 

190.9 

6.945 

4.85 

O?  (very  crowded  field,  small  PM) 
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WDS 

(2000) 

HIPl 

Discoverer 

Designation 

N 

0 

(deg) 

p 

(") 

Am 

(mag) 

Status,  remarks 

19234+2034 

95309 

RA068AB 

1 

350.9 

5.13 

4.7 

P  (2MASS  position) 

19287+3437 

95769 

RA0  86AB 

1 

348.2 

9.916 

5.75 

0  (2MASS  position,  high  crowding) 

19359+5659 

96395 

RA0  87AB 

1 

250.7 

10.08 

5.3 

P  (2MASS  position  &  color) 

19464+3344 

97222 

RA069FK 

1 

318.3 

16.402 

6.21 

0  (2MASS  position) 

20086+8507 

99232 

RAO  21  AC 

2 

65.3 

1.32 

5.5 

P?  (small  separation,  low  crowding) 

20124-1237 

99572 

RAO  88  CD 

1 

271.1 

17.823 

5.14 

0  (not  found  in  2MASS) 

20169+5017 

99965 

RAO70FL 

1 

264.2 

15.597 

4.66 

0  (2MASS  position,  high  crowding) 

20312+5653 

101234 

RA022AB 

1 

195.8 

0.175 

0.0 

P  (acceleration  binary,  conf.  with  P200) 

20333+3323 

101430 

RA071  AE 

1 

224.7 

12.279 

6.18 

0  (2MASS  position,  high  crowding) 

20408+1956 

102040 

RAO  23  Ca,Cb 

1 

59.4 

0.266 

0.0 

P  (close  and  bright) 

20408+1956 

102040 

RAO  23  CD 

1 

303.0 

20.326 

3.84 

0  (not  found  in  2MASS,  high  crowding) 

20408+1956 

102040 

RAO  23  CE 

1 

327.2 

29.350 

4.37 

0  (not  found  in  2MASS,  high  crowding) 

20577+2624 

103455 

RA024AB 

2 

102.8 

0.610 

4.9 

P  (acceleration  binary,  conf.  with  P200) 

20599+4016 

103641 

RA072DE 

1 

117.6 

19.752 

5.73 

0  (2MASS  position,  crowding),  D  =  HIP  103052 

21102+2045 

104514 

RA0  25AB 

2 

209.1 

3.39 

4.6 

P  (2MASS  image,  conf.  with  P200) 

21585+0347 

108473 

RA073AB 

1 

89.3 

12.39 

5.6 

P  (2MASS  position  &  color) 

22094+3508 

109361 

RA026AB 

1 

351.0 

0.368 

2.5 

P  (aceleration  &  SBl) 

22204+4625 

110291 

RA074AB 

1 

122.9 

14.980 

4.52 

0  (2MASS  position,  crowding) 

22204+4625 

110291 

RAO  74  AC 

1 

232.1 

30.572 

3.82 

0  (2MASS  position,  crowding) 

22240+0612 

110574 

RA027Aa,Ab 

1 

327.0 

0.092 

0.0 

P  (acceleration  binary,  conf.  with  P200) 

22246+3926 

110626 

RA028AB 

1 

315.8 

4.37 

6.0 

P  (2MASS  image) 

22311+4509 

111148 

RA029BC 

1 

215.7 

3.9108 

6.00 

0  (not  found  in  2MASS,  high  crowding) 

22311+4509 

111148 

RA029BD 

1 

264.4 

12.276 

3.59 

0  (2MASS  position,  crowding) 

22524+0950 

112935 

RAO  30  Da, Db 

4 

319.3 

0.242 

0.2 

P  (close,  fixed  for  1  yr) 

23108+4531 

114456 

RAO  89  CD 

1 

22.3 

20.587 

2.63 

0  (2MASS  position) 

23108+4531 

114456 

RAO  89  CE 

1 

290.6 

15.330 

3.03 

0  (2MASS  position) 

23258+4521 

115655 

RA0  75BC 

1 

65.4 

13.195 

5.79 

0?(2MASS  color) 

23258+4521 

115655 

RAO  75  CD 

1 

279.8 

0.612 

1.01 

0  (unrelated  to  B) 

23419-0559 

116906 

RA0  31Ba,Bb 

3 

69.6 

0.519 

0.1 

P  (close  &  bright,  fixed  for  1  yr) 

23588+3156 

118213 

RA076AB 

2 

87.8 

4.82 

6.3 

P  (2MASS  image,  conf.  with  P200) 

23588+3345 

118225 

RA077AB 

2 

172.3 

4.91 

7.0 

P  (2MASS  image,  fixed  for  1  yr) 

separation,  color,  crowdedness  of  the  field,  and  proper  motion 
of  the  main  target.  We  do  not  compute  the  formal  probability 
of  physical  association  based  on  relative  astrometry  because  the 
motion  can  be  distorted  by  subsystems,  while  chance  projections 
with  small  and  similar  motions  do  happen  sometimes  (see 
e.g.,  Tokovinin  &  Lepine  2012).  Other  factors  involved  in  the 
classification  are  difficult  to  quantify  in  terms  of  probability. 
Most  faint  companions  in  Table  5  wider  that  10''  are  optical 
with  a  large  confidence;  only  three  of  them  have  an  uncertain 
status.  In  some  instances,  the  companions  are  not  found  in  the 
catalog,  but  the  2MASS  images  in  the  Ks  band  reveal  their 
presence  as  blends.  This  is  taken  as  confirmation,  because 
the  companion  position  matches  (at  least  qualitatively)  and 
because  the  detection  of  the  blend  in  2MASS  indicates  that 
the  companion  is  brighter  in  Ks  than  in  /',  indicating  it  is  a 
cooler  object  below  the  detection  threshold  in  the  catalog. 

Notes  in  the  Appendix  give  additional  information  on  the 
multiple  systems  appearing  in  Table  5,  the  status  of  the  measured 
binary  pairs  (“P”  for  true  physical  binaries,  “O”  for  chance 
optical  alignments  of  unrelated  stars),  and  the  reasons  for 
this  classification.  Follow-up  observations  in  September  2013 
with  the  PALM  3000  AO  system  at  the  5  m  Hale  telescope 
at  Palomar  Observatory  (Dekany  et  al.  2013)  are  mentioned 
where  relevant,  while  their  full  results  will  be  published  in  a 
forthcoming  paper  (L.  C.  Roberts  et  al.,  in  preparation).  Two 
of  the  systems  with  newly  resolved  binaries,  HIP  2292B  (= 
HIP  2350)  and  HIP  121 89B  (=  HIP  12184),  are  also  exoplanet 
hosts;  a  complete  discussion  of  their  properties  are  included  in 
Roberts  et  al.  (2014).  Two  RAO  pairs  were  also  measured  by 


speckle  interferometry  at  the  SOAR  4.1  m  telescope  (Tokovinin 
et  al.  2014). 

5.2.  Binarity  of  Secondary  Components 

Overall,  there  are  212  secondary  components  with  separations 
above  10"  that  were  observed  with  Robo-AO.  Figure  6  illustrates 
different  hierarchies  found  in  this  subsample.  The  wide  binary 
is  the  root  of  the  hierarchy  (level  1);  in  112  cases  there  are  no 
subsystems.  Another  100  wide  binaries  contain  inner  hierarchies 
of  level  11  (subsystem  in  the  primary  component),  level  12 
(subsystem  in  the  secondary)  or  both  (that  is,  a  2-^2  quadruple). 
Here  we  silently  ignore  several  subsystems  (levels  1 1 1  etc.),  and 
restrict  the  discussion  to  the  hierarchical  levels  1,11,  and  12.  As 
some  subsystems  are  most  likely  not  yet  discovered,  it  is  safe  to 
say  that  at  least  half  of  those  wide  binaries  are  in  fact  triple  or 
higher-order  multiples. 

Orbital  periods  of  wide  binaries  are  estimated  roughly  from 
their  projected  separation  and  denoted  as  P*  to  distinguish 
them  from  orbital  solutions.  Masses  of  binary  components  and 
the  mass  ratio  {q  =  M2IM1)  are  estimated  from  their  ab¬ 
solute  magnitudes,  assuming  that  the  stars  obey  standard  re¬ 
lations  for  the  main  sequence.  The  detection  limits  of  imag¬ 
ing  are  converted  from  the  observed  parameters  (p.  Am)  into 
binary  parameters  (P*,^),  and  combined  with  the  detection 
limits  from  spectroscopy  and  other  techniques  (more  details  in 
Tokovinin  2014). 

There  are  83  primary  (level  11)  and  39  secondary  (level  12) 
subsystems  in  the  surveyed  objects  (Figure  7);  19  new  secondary 
subsystems  were  found,  doubling  their  previously  observed 
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Figure  6.  Observed  numbers  of  different  hierarchies  in  212  wide  binaries  surveyed  by  Robo-AO. 


Separation,  arcsec 

0.001  0.010  0.100  1.000  10.000  100.000 


Period,  d 


Figure  7.  Distribution  of  39  secondary  subsystems  (upper  panel)  and  83  primary 
subsystems  (lower  panel)  of  wide  binaries  in  the  {P,q)  plane.  The  contours 
indicate  average  detection  probability.  The  angular  separation  on  the  upper  axis 
corresponds  to  the  distance  of  50  pc. 


fraction.  At  first  glance,  the  primary  subsystems  still  dominate, 
but  remember  that  most  of  the  primary  targets  (bright  F  and 
G  dwarfs)  were  surveyed  in  radial  velocity  (RV),  while  only  a 
few  brighter  secondaries  have  such  RV  coverage.  The  primaries 
also  benefit  from  the  Hipparcos  discoveries  of  acceleration 


Table  6 

Frequency  of  Levels  11  and  12  with  3.5  ^  Xm  ^5 


Parameter 

Primary 

Secondary 

Total  number 

25 

20 

Raw  frequency 

0.118  ±0.023 

0.094  ±  0.021 

Average  detection 

0.88 

0.58 

Frequency  =  0) 

0.153  ±  0.003 

0.164  ±0.037 

Frequency  {^  =  \) 

0.130  ±0.026 

0.127  ±0.028 

subsystems  (Makarov  &  Kaplan  2005),  which  explains  deeper 
detections  at  P  10^  d.  At  P  10^  d  the  situation  is  reversed  in 
favor  of  secondaries,  owing  to  their  imagery  with  Robo-AO.  At 
still  longer  periods  (separation  the  companion  census  for 
primaries  is  complete  because  of  2MASS  (Tokovinin  2011)  and 
CPM  (Tokovinin  &  Lepine  2012)  data,  while  for  the  secondaries 
the  detections  are  restricted  by  the  size  of  the  Robo-AO  field. 
However,  this  makes  little  difference  because  extremely  wide 
subsystems  are  dynamically  unstable. 

It  is  convenient  to  express  orbital  periods  by  v,  the  logarithm 
of  the  period  in  days.  The  criterion  of  dynamical  stability  by 
Mardling  &  Aarseth  (2001)  requires  that  the  ratio  of  outer  to 
inner  periods  in  triple  systems  be  greater  than  4.7,  depending 
on  the  eccentricity  in  the  outer  orbit.  The  dynamical  truncation 
function  P(Av)  (where  Ax  =  Xout  —  -^in  is  the  logarithm  of  the 
period  ratio)  is  the  probability  of  a  given  triple  system  being 
dynamically  stable;  it  is  modeled  here  as  zero  for  Ax  <  0.7,  one 
for  Ax  >  1.7,  and  linear  in  between.  Considering  the  large  range 
in  periods,  the  exact  form  of  this  function  has  little  influence  on 
the  results. 

In  the  subsample  of  212  wide  binaries,  the  average  dynamical 
truncation  is  0.95  at  inner  period  vin  =  5  and  0.5  at  Xin  =  6, 
explaining  the  lack  of  wide  subsystems  at  levels  1 1  and  12.  There 
is  no  need  to  look  for  binary  secondaries  with  separations  >  10'' 
because  those  are  intrinsically  rare.  The  limitation  imposed 
by  the  Robo-AO  field  of  view  is  therefore  not  important  for 
this  work. 

We  select  for  analysis  the  period  range  3.5  ^  Xin  ^  5  where 
the  dynamical  truncation  is  not  important  and  the  detection 
of  secondary  subsystems  with  Robo-AO  is  relatively  complete. 
This  range  corresponds  to  one  decade  in  separation.  We  compare 
the  frequency  of  primary  and  secondary  subsystems  in  this 
range  of  periods  among  selected  212  wide  binaries  in  Table  6. 
The  number  of  missed  companions  obviously  depends  on 
the  distribution  of  the  mass  ratio,  modeled  as  a  power  law 
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f{q)  a  It  is  generally  accepted  that  the  power  index  ^ 
for  solar-type  binaries  is  close  to  zero  (uniform  distribution) 
or  slightly  negative  (Duchene  &  Kraus  2013),  while  low-mass 
binaries  tend  to  have  more  equal  components  and  ^0  ~  1 . 

The  probability  of  detecting  a  subsystem  averaged  over  period 
(in  the  selected  interval)  and  over  the  mass  ratio  q  can  be  used 
to  correct  the  raw  companion  frequency,  if  ^0  =  0.  However,  the 
data  collected  here  indicate  that  in  the  secondary  subsystems 
^  ^  requiring  a  smaller  correction  for  missed  binaries. 
Therefore,  the  estimated  frequency  of  secondary  subsystems 
with  3.5  ^  Vin  ^  5  is  0.13  zb  0.03,  assuming  conservatively 
P  =  1.  It  is  only  slightly  less  than  the  frequency  of  primary 
(level  11)  subsystems,  which  is  less  sensitive  to  the  assumed  ^ 
and  is  somewhere  between  0.13  and  0.15  (Table  6). 

The  frequency  of  subsystems  found  here  is  close  to  the 
frequency  of  solar-type  binaries  in  the  same  period  range.  For 
the  log-normal  period  distribution  and  the  binary  fraction  of  0.46 
derived  by  Raghavan  et  al.  (2010),  the  fraction  of  companions 
with  3.5  ^  Vin  ^  5  is  0.112  (a  similar  estimate  is  obtained 
for  the  full  FG-67  sample).  The  frequency  of  subsystems  in  the 
components  of  wide  binaries  may  be  enhanced  in  comparison 
to  the  single  stars,  but  only  slightly. 

The  actual  frequency  of  subsystems  in  secondaries  may  be 
even  higher  because  of  two  additional  factors.  First,  the  methods 
used  to  identify  wide  secondaries  introduce  some  selection 
against  close  binaries  (Tokovinin  &  Lepine  2012).  Second, 
detection  limits  here  may  be  over-estimated  (see  Section  3.4). 
When  the  assumed  detection  limits  are  deeper  than  in  reality, 
the  correction  of  the  raw  frequency  for  incomplete  detection 
becomes  smaller. 

Quite  surprisingly,  the  subsystems  in  primary  and  secondary 
are  often  found  at  the  same  time.  Among  the  39  secondary 
subsystems,  22  (more  than  half!)  belong  to  2-^2  quadruples 
(Figure  6).  Considering  that  many  secondary  subsystems  with 
short  periods  are  still  missed,  the  actual  fraction  of  2-^2  quadru¬ 
ples  among  212  wide  binaries  should  be  more  than  10%.  This 
surprising  correlation  likely  has  a  root  in  stellar  system  forma¬ 
tion  processes  and  requires  further  study. 

5.3.  Frequency  of  Tertiary  Components 

In  the  second  part  of  our  survey  we  looked  for  distant 
companions  to  close  binaries.  In  other  words,  the  problem  is 
“reversed”  and  we  explore  the  hierarchy  from  inside-out.  Similar 
studies  were  made  by  Tokovinin  et  al.  (2006),  Allen  (2007), 
Rucinski  et  al.  (2007),  and  others. 

The  second  survey  selects  main  targets  that  are  close  binaries 
with  P  <  100  yr  (x  <  4.56).  Binaries  with  yet  unknown 
periods,  such  as  acceleration  binaries  from  Hipparcos  (Makarov 
&  Kaplan  2005)  and/or  spectroscopic  binaries  from  Nordstrom 
et  al.  (2004),  are  included.  We  resolved  three  acceleration 
binaries  (HIP  6653,  101234,  103455)  with  Robo-AO,  while 
more  were  resolved  by  a  targeted  campaign  with  the  NICI 
instrument  (Tokovinin  et  al.  2012,  2013).  However,  as  both 
acceleration  and  RV  techniques  involve  a  non-negligible  false 
alarm  probability,  some  of  those  presumed  close  binaries  are  in 
fact  single  stars. 

For  the  purpose  of  statistical  analysis,  we  restrict  the  sample 
to  241  inner  binaries  with  known  periods  P  <  lOOyr.  In  the 
cases  where  even  closer  inner  subsystems  are  present  (e.g., 
the  A-component  of  a  visual  binary  is  a  spectroscopic  pair),  we 
select  the  inner-most  (closest)  pair  and  proceed  outward.  The 
71  close  binaries  have  additional  outer  (tertiary)  companions; 
some  of  the  tertiary  components  are  themselves  close  binaries 


Separation,  arasea 

0.001  0.010  0.100  1.000  10.000  100.000 


Figure  8.  Periods  and  mass  ratios  of  71  tertiary  components  to  241  binaries  with 
P*  <  lOOyr  observed  with  Robo-AO.  The  contours  indicate  average  detection 
probability.  The  full  line  over-plots  the  average  dynamical  truncation  (not  to  be 
confused  with  q). 


(in  2+2  quadruples),  and  sometimes  there  are  even  more  distant 
companions  (3-fl  quadruples).  The  raw  frequency  of  tertiaries 
is  29%.  When  we  restrict  the  sample  to  shorter  inner  periods, 
the  fraction  of  triples  increases,  but  the  sample  becomes  smaller 
and  the  statistical  errors  increase. 

Figure  8  presents  tertiary  components  in  the  (Pout.  Qi)  plane, 
where  q\  =  Afs/Af  i  is  the  ratio  of  the  masses  between  tertiary 
and  primary  components.  The  iso-detection  curves  go  deeper 
than  for  the  general  sample  (compare  with  Figure  7,  lower  panel) 
owing  to  the  Robo-AO  imaging.  Many  of  tertiary  components 
are  found  in  the  2MASS  catalog,  but  their  colors  alone  were 
not  sufficient  to  establish  that  they  are  physical,  especially  in 
crowded  fields.  The  second-epoch  imaging  with  Robo-AO  helps 
here.  A  few  uncertain  cases  (companions  detected  by  Robo-AO 
in  crowded  fields  without  2nd  epoch  data)  are  not  yet  accepted 
as  real. 

The  average  dynamical  truncation  curve  P(xout)  is  over¬ 
plotted.  It  explains  why  tertiaries  with  short  periods  Xout  ^  3  are 
rare.  Such  tertiaries  would  be  easily  detectable  by  spectroscopy 
and  their  scarcity  is  genuine.  As  in  Section  5.2,  we  evaluate  the 
frequency  of  outer  systems  in  the  outer  period  range  of  1.5  dex, 
6  ^  -^out  ^  2.5 — the  range  free  from  dynamical  truncation  and 
well  covered  by  the  detection  techniques.  It  is  safe  to  assume 
f  =  0  for  wide  components,  in  which  case  the  average  detection 
probability  is  0.88.  The  raw  companion  frequency  of  34/241  = 
0.141  zb  0.024  is  corrected  to  0.161  zb  0.028.  Owing  to  the 
deep  detection  limit,  there  is  little  difference  between  the  raw 
and  detection-corrected  estimates.  The  frequency  is  significantly 
higher  than  the  companion  frequency  to  main  targets  in  the 
same  period  range,  0.090  (again  using  the  data  of  Raghavan 
et  al.  2010). 

It  is  premature  to  speculate  about  the  meaning  of  this  finding. 
In  Figure  8  we  note  the  “pile-up”  of  tertiary  companions 
at  Vout  2  (separations  on  the  order  of  10''  and  projected 
separations  on  the  order  of  500  AU).  If  this  feature  is  not  a 
statistical  fluctuation,  it  might  imply  some  characteristic  scale 
in  the  formation  of  solar-type  multiple  stars. 

The  overall  raw  frequency  of  tertiary  components  (with  all 
periods)  is  only  29%.  If  we  assume  that  the  intrinsic  distribution 
of  Vout  is  Gaussian  (Raghavan  et  al.  2010)  and  that  their  mass 
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Table  7 

New  Orbital  Elements 


WDS 

HIPl 

Figure 

P 

(yr) 

a 

(") 

i 

(deg) 

Q 

(deg) 

To 

(yr) 

e 

CO 

(deg) 

Gr 

Ref. 

01158+0947 

5898 

9a 

162. 

0.607 

101.7 

151.2 

2101. 

0.447 

262.8 

4 

a 

±13. 

±0.035 

±1.0 

±2.1 

±13. 

±0.054 

±6.6 

01417-1119 

7916 

9b 

601. 

2.42 

82.6 

267.98 

2015.23 

0.901 

348.3 

4 

b 

±107. 

±0.19 

±2.2 

±0.93 

±0.87 

±0.030 

±9.7 

02371-1112 

12204 

9c 

351. 

1.01 

64.4 

277. 

1902. 

0.508 

309. 

5 

c 

±271. 

±0.34 

±3.4 

±12. 

±15. 

±0.044 

±43. 

03140+0044 

15058 

9d 

753. 

1.32 

132.5 

123.7 

1883. 

0.281 

254. 

4 

d 

±84. 

±0.13 

±2.8 

±7.8 

±27. 

±0.020 

±30. 

03413+4554 

17217 

10a 

46.2 

t)21 

90.3 

269.9 

1974.8 

0.79 

50.1 

4 

±1.1 

±0.13 

±3.6 

±4.5 

±1.7 

±0.76 

±9.0 

03496-0220 

17895 

10b 

54. 

0.49 

114.5 

234. 

1992.8 

0.750 

105. 

4 

±29. 

±0.12 

±1.4 

±23. 

±3.4 

±0.012 

±12. 

04163+0710 

19911 

10c 

12.38 

0.2016 

69.99 

144.8 

2004.40 

0.1677 

65.6 

3 

±0.03 

±0.0066 

±0.75 

±1.5 

±0.04 

±0.0043 

±1.5 

04422+3731 

21878 

lOd 

564. 

1.51 

141.2 

209.1 

2024.9 

0.4531 

262. 

4 

e 

±119. 

±0.22 

±6.0 

±7.6 

±8.0 

±0.0071 

±11. 

07043-0303 

34110 

11a 

46.88 

0.273 

99.2 

99.4 

1968.59 

0.635 

33.1 

3 

f 

±0.50 

±0.014 

±2.2 

±1.1 

±0.76 

±0.037 

±5.9 

07092+1903 

34524 

11b 

27.60 

0.179 

103.3 

355.6 

2010.81 

0.94 

244.5 

4 

±1.1 

±0.013 

±1.9 

±1.7 

±0.25 

* 

±2.3 

08122+1739 

40167 

11c 

19.05 

0.3728 

149.2 

262. 

1983.05 

0.088 

41. 

4 

g 

±0.46 

±0.0068 

±7.2 

±12. 

±0.32 

±0.013 

±23. 

10387+0544 

52097 

lid 

488.0 

1.1962 

66. 

159.89 

1748. 

0.5967 

339.20 

4 

±3.6 

±0.0051 

±22. 

±0.36 

±14. 

±0.0018 

±0.29 

11293+3025 

56054 

12a 

254.9 

1.18 

117. 

246. 

2075. 

0.65 

92. 

5 

±7.9 

±0.12 

±98. 

±20. 

±61. 

±0.12 

±48. 

13577+5200 

68193 

12b 

205. 

0.857 

139.5 

276.8 

1911.6 

0.717 

330. 

3 

h 

±20. 

±0.028 

±9.2 

±5.5 

±2.3 

±0.025 

±11. 

15413+5959 

76837 

12c 

434. 

1.013 

74.4 

211.83 

1959.42 

0.828 

21.4 

4 

i 

±30. 

±0.036 

±1.2 

±0.74 

±0.95 

±0.016 

±2.2 

18154+5720 

89455 

12d 

82. 

0.40 

62. 

259. 

1997.8 

0.36 

313. 

4 

±65. 

±0.21 

±14. 

±19. 

±6.6 

±0.36 

±51. 

20514-0538 

102945 

13a 

199.8 

0.844 

65.13 

173.36 

1896.99 

0.526 

47.9 

2 

j 

±3.0 

±0.015 

±0.70 

±0.86 

±0.81 

±0.014 

±2.4 

21041+0300 

103987 

13b 

21.5 

0.231 

65.4 

102.3 

1985.0 

0.24 

25. 

3 

±6.6 

±0.041 

±6.2 

±6.6 

±6.9 

±0.21 

±18. 

21243+3740 

105676 

13c 

17.8 

0.167 

122.1 

247.3 

2014.6 

0.31 

285.6 

3 

±1.9 

±0.015 

±4.8 

±7.8 

±1.3 

±0.10 

±7.4 

22161-0705 

109951 

13d 

44. 

0.30 

36. 

160. 

1989.0 

0.55 

126. 

4 

k 

±24. 

±0.16 

±101. 

±62. 

±5.3 

±0.59 

±104. 

References,  (a)  Docobo  &  Ling  (2010);  (b)  Rica  (2013);  (c)  Heintz  (1998);  (d)  Heintz  (1963);  (e)  Mason  et  al.  (2004a);  (f)  Docobo  &  Ling  (2009);  (g)  Heintz  (1996b); 
(h)  Heintz  (2001);  (i)  Heintz  (1975);  (j)  Heintz  (1998);  (k)  Cvetkovic  (2011) 


ratio  qi  is  distributed  uniformly,  the  companion  frequency  of  this 
underlying  distribution  (before  applying  dynamical  truncation 
and  detection  filters)  is  determined  to  be  0.54;  i.e.,  slightly 
but  not  dramatically  enhanced  in  comparison  with  0.46  for  all 
solar-type  dwarfs.  Most  binaries  with  P  <  lOOyr  do  not  have 
outer  tertiary  companions;  a  similar  conclusion  was  reached  by 
Tokovinin  et  al.  (2006)  for  spectroscopic  binaries  with  periods 
of  "^10  d  and  longer.  On  the  other  hand,  the  closest  binaries, 
with  P  <  3  d,  are  found  almost  exclusively  in  triples. 

5.4.  New  and  Updated  Orbits 

For  purposes  of  orbit  determination,  Robo-AO  provides  a 
unique  collection  of  data.  While  there  are  certainly  exceptions 
(Tokovinin  2012;  Horch  et  al.  2012),  the  majority  of  interfero¬ 
metric  data,  whether  with  a  filled  or  dilute  aperture,  have  modest 
differential  magnitude  limits  of  about  Am  =  3  mag.  For  many 
pairs  of  larger  magnitude  difference,  the  Robo-AO  observation 
represent  the  first  observation  by  a  high-resolution  technique. 


and  given  the  decline  in  classical  micrometry,  the  first  obser¬ 
vation  of  any  kind  in  many  decades.  Given  these  factors,  the 
number  of  first  orbits  as  well  as  those  with  significant  changes 
from  the  previous  determinations  is  not  surprising. 

New  orbits  are  presented  in  Table  7.  In  this  table,  the  pairs  are 
identified  by  their  WDS  and  HIPl  numbers,  followed  by  their 
orbital  elements,  giving  the  period  P  in  years,  the  semimajor  axis 
a  in  arcseconds,  the  inclination  i  and  longitude  of  the  node  Q, 
both  in  degrees,  the  epoch  of  the  most  recent  periatron  passage 
To  in  Besselian  years,  the  eccentricity  e  and  the  longitude 
of  periastron  co  in  degrees.  Following  this  is  the  orbit  grade 
(see  Hartkopf  &  Mason  2013b  for  details)  as  an  evaluation  of 
the  orbit.  For  those  pairs  with  previous  orbit  solutions  we  provide 
a  reference  to  the  previous  “best”  orbit.  Figures  9-13  show  the 
new  orbital  solutions,  plotted  with  all  published  data  in  the  WDS 
database.  In  each  of  these  figures,  micrometric  observations 
are  indicated  by  plus  signs,  high  resolution  measures  by  filled 
circles,  Robo-AO  observations  as  filled  stars,  and  Hipparcos 
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Figure  9.  New  orbits  for  the  systems  listed  in  Table  7,  together  with  all  published  data  in  the  WDS  database.  See  Section  5.4  for  a  description  of  symbols  used  in  this 
and  the  following  figures.  In  all  of  these  figures,  panel  (a)  is  the  upper  left  corner,  (b)  the  upper  right,  (c)  the  lower  left,  and  (d)  the  lower  right,  (a)  HIP  5898;  (b)  HIP 
7916;  (c)  HIP  12204;  (d)  HIP  15058. 


observations  as  filled  diamonds.  An  O— C  line  connects  each 
measure  to  its  predicted  position  along  the  new  orbit  (shown  as 
a  thick  solid  line).  The  shaded  circle  centered  on  the  primary 
indicates  the  resolution  limit  of  the  P60.  A  dot-dashed  line 
indicates  the  line  of  nodes,  and  a  curved  arrow  in  the  lower 
right  comer  of  each  figure  indicates  the  direction  of  orbital 
motion.  Previous  published  orbits  are  shown  as  dashed  ellipses; 
references  to  each  of  the  published  orbits  are  given  in  the  notes 
to  Table  7. 

Orbital  elements  were  determined  with  orbgridlO,  the  latest 
version  of  the  venerable  orbit  reduction  package  of  Hartkopf 
et  al.  (1989),  which  utilizes  a  three  dimensional  (P,  T,  and  e) 
grid  search  of  variable  grid  step  sizes  from  three  initial  estimates. 


It  has  been  modified  to  include  weighting  methodologies  deter¬ 
mined  in  the  production  of  the  5*  Catalog  of  Orbital  Elements 
of  Visual  Binary  Stars  (Hartkopf  et  al.  2001).  Some  orbits  use 
new  speckle  data  from  the  SOAR  telescope  (Tokovinin  et  al. 
2014). 

The  rest  of  this  section  gives  notes  on  some  pairs  with  new 
orbits. 

HIP  12204  (=WDS  J02371-1112  =  HU  1216AB)  has  only 
a  preliminary  orbit  of  grade  5,  hence  the  formal  errors  of  the 
elements  are  very  large. 

HIP  17895  (=WDS  J03496-0220  =  YR  23):  This  is  the  first 
calculated  orbit  of  this  pair  announced  in  Horch  et  al.  (2002b) 
based  on  observations  with  the  WIYN  telescope.  The  2004 
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Figure  10.  Same  as  Figure  9.  (a)  HIP  17217;  (b)  HIP  17895;  (c)  HIP  19911;  (d)  HIP  21878. 


measure,  obtained  with  the  26  inch  telescope  in  Washington, 
D.C.  (Mason  et  al.  2006),  is  off  a  bit,  which  is  not  surprising 
given  the  aperture. 

HIP  34524  (=WDS  J07092+1903  =  CHR  216)  was  resolved 
by  lunar  occultations  by  Africano  et  al.  (1975)  and  has  several 
speckle  measurements.  Robo-AO  finds  it  at  nearly  the  same  po¬ 
sition  as  Hipparcos  (the  quadrant  of  our  measure  is  established 
by  the  automatic  data  reduction  co-added  image),  so  the  binary 
has  made  nearly  one  full  revolution  since  1991.25.  It  was  not 
resolved  at  SOAR  in  2009-2011  (p  <  0'.'03),  but  measured  in 
2012.9,  2013.05,  and  2014.06  (these  data  are  used  in  the  orbit 
calculation).  The  orbit  is  very  eccentric,  but  measurements  do 
not  constrain  the  eccentricity,  so  it  was  fixed  to  0.94 — a  provi¬ 
sional  value  that  gives  the  expected  mass  sum  of  1.8  AIq-  The 


orbit  should  be  followed  by  spectroscopy.  Unfortunately,  the 
periastron  passage  in  2010  was  missed,  we  have  to  wait  for  the 
next  one  in  2038. 

HIP  40167  (=WDS  J08122-F1739  =  HUT  lCa,Cb):  This  is 
the  first  calculated  orbit  of  this  pair.  The  existence  of  the  Ca,Cb 
pair  was  inferred  and  the  first  astrometric  orbit  calculated  by 
Heintz  (1996b). 

HIP  103987  (=WDS  121041-^0300  =  WSI 6):  This  is  the  first 
calculated  orbit  of  this  pair  first  split  by  Mason  et  al.  (2001a) 
with  the  Struve  2.1m  telescope  in  an  investigation  of  Hipparcos 
acceleration  (G-type  solution)  pairs. 

HIP  105676  (=WDS  121243-^3740  =  WSI  7):  Like  HIP 
103987  above,  this  is  the  first  calculated  orbit  of  this  pair  also 
first  resolved  with  the  same  telescope  on  the  same  project. 
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Figure  11.  Same  as  Figure  9.  (a)  HIP  341 10;  (b)  HIP  34524;  (c)  HIP  40167;  (d)  HIP  52097. 


6.  SUMMARY 

A  survey  of  high-resolution  imaging  with  Robo-AO  has 
advanced  our  knowledge  of  the  multiplicity  of  solar-type  stars 
in  several  ways. 

The  binarity  of  the  wide  secondary  components  turned  out  to 
be  comparable  to  the  frequency  of  subsystems  in  the  primary 
components.  This  overturns  the  established  paradigm  about 
multiple-star  architecture,  where  subsystems  were  thought  to 
be  preferentially  associated  with  the  primary,  most  massive  star 
in  the  system.  This  traditional  view  can  be  traced  to  the  previous 
technological  restrictions  allowing  to  study  only  bright  stars  by 
spectroscopy  or  high  angular  resolution,  and  matched  the  idea 
that  multiple  stars  experience  chaotic  dynamical  interactions 
where  the  least  massive  components  are  ejected  onto  distant 
orbits  or  leave  the  system  altogether.  This  new  paradigm  of 
dynamical  evolution  of  multiple  star  systems,  where  wide 


secondary  components  retain  companions  at  the  same  frequency 
as  the  primary,  will  require  more  exploration  to  determine  how 
such  systems  develop. 

Quite  unexpectedly,  we  discovered  that  subsystems  in  the 
primary  and  secondary  components  of  wide  binaries  are  of¬ 
ten  present  simultaneously,  i.e.,  these  systems  are  in  fact  2-\-2 
quadruples.  This  new  result  means  that  2-\-2  quadruples  are  rel¬ 
atively  frequent.  Note  that  in  the  25  pc  sample  of  Raghavan 
et  al.  (2010),  9  out  of  11  quadruples  have  the  2-^2  architec¬ 
ture  and  only  2  are  of  the  3-^1  “planetary”  configuration;  this 
result,  based  in  small-number  statistics,  is  now  supported  by 
the  much  larger  FG-67  sample  studied  with  Robo-AO.  The 
high  frequency  of  2-^2  quadruples  challenges  the  paradigm 
of  chaotic  A-body  dynamics  (which  can  produce  such  ar¬ 
chitecture  only  exceptionally)  and  calls  for  exploration  of 
other  formation  scenarios  for  multiple  stars  (e.g.,  Whitworth 
2001). 
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Figure  12.  Same  as  Figure  9.  (a)  HIP  56054;  (b)  HIP  68193;  (c)  HIP  76837;  (d)  HIP  89455. 


The  high  dynamic  range  of  the  Robo-AO  imaging  enabled 
discovery  of  17  certain  and  2  possible  tertiary  components 
to  close  binaries,  converting  them  into  triple  or  higher-order 
systems.  Statistics  of  such  triples  reveal  some  interesting  details, 
such  as  the  prevalence  of  outer  separations  on  the  order  of 
500  AU.  This  result  is  not  yet  formally  significant,  but  it 
highlights  the  potential  of  large  samples.  Historically  small 
binary  samples  allowed  only  the  first-order  description  of  binary 
statistics  by  smooth  functions  with  a  few  parameters  (Duchene 
&  Kraus  2013).  Modern  large  and  homogeneous  samples  begin 
to  reveal  new  details  in  these  distributions,  such  as  twin  binaries 
with  identical  masses  (Lucy  2006),  bimodal  period  distribution 
in  Hyades  (Griffin  2012),  and,  possibly,  preferred  separations  of 
tertiary  components  seen  here.  These  findings  will  advance  our 
understanding  of  multiple-star  formation. 


The  large  data  set  that  resulted  from  this  survey  enabled  us  to 
compute  first  visual  orbits  for  9  pairs  and  to  update  (sometimes 
dramatically)  the  existing  orbits  for  11  more.  Accumulation 
and  improvement  of  the  visual-orbit  data  will  lead  to  a  better 
knowledge  of  stellar  masses  and  to  the  statistical  analysis  of 
orbital  elements.  It  will  also  enable  dynamical  study  of  visual 
binaries  presenting  special  interest. 
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Figure  13.  Same  as  Figure  9.  (a)  HIP  102945;  (b)  HIP  103987;  (c)  HIP  105676;  (d)  HIP  109951. 
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APPENDIX 

NOTES  ON  NEW  BINARY  STAR  SYSTEMS 

HIP  2292B  =  HIP  2350  =  RAO  IBC  is  resolved  at  (T.'S,  with 
several  Robo-AO  measures  and  a  confirmation  with  the  Palomar 
200  inch  (5  m)  telescope  (P200).  It  remained  fixed  during  1  yr, 
despite  the  proper  motion  (PM)  of  0'.'25  yr“^ .  This  star  also  hosts 
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an  exoplanet  system  with  P  =  3.44  d.  The  common-proper- 
motion  (CPM)  component  A  is  at  839''  from  B.  We  discuss  this 
multiple  system  in  a  separate  paper. 

HIP  2717B  =  BD+ei  119(7  =  8.98,  F8V,  X-ray  source)  has 
two  new  companions  at  ~3",  designated  as  RAO  2BC  and  RAO 
2BD.  We  see  the  brighter  companion  C  in  the  2MASS  image 
at  a  similar  position,  although  it  is  not  listed  as  a  separate  point 
source  in  2MASS.  On  the  other  hand,  the  fainter  star  D  is  likely 
optical  because  the  field  is  crowded  and  because  a  trapezium¬ 
like  configuration  of  BC  and  BD  would  be  dynamically  un¬ 
stable.  The  CPM  nature  of  AB  (Tokovinin  &  Lepine  2012) 
is  confirmed  by  the  data  from  2MASS  and  UCAC  (Zacharias 
etal.  2013). 

HIP  3214B  (V  =  14.87)  is  resolved  at  1'.'6.  The  new  pair 
is  considered  physical  because  the  companion  is  bright  and 
because  it  remained  fixed  during  1  yr.  The  2MASS  image  is 
elongated  in  the  E-W  direction. 

HIP  3540  AB  is  LDS  836  at  55'.'4  separation.  We  resolved  B 
into  a  new  1'.'5  pair  RAO  4BC,  confirmed  with  the  P200.  The 
component  A  is  a  spectroscopic  binary  (SB)  with  P  =  11.4  d 
(Griffin  2002).  The  system  is  thus  a  2-^2  quadruple. 

HIP  3795  is  another  new  quadruple  system.  The  AB  pair 
BU  495  has  an  orbit  with  P  =  413  yr.  The  CPM  component 
C  (7  =  10.83'”)  at  152"  from  AB  turns  out  to  be  a  new  O'.' 67 
pair  RAO  5Ca,Cb.  We  consider  Ca,Cb  physical  because  the 
companion  is  bright,  the  separation  is  small,  and  the  field  is  not 
crowded. 

HIP  401 6B  is  accompanied  by  the  1'.'7  pair  CD,  at  a  distance 
of  17"  from  B  (AB  =  DAL  11  at  41"  is  physical).  The  field  is 
crowded,  so  we  consider  the  faint  stars  C  and  D  to  be  unrelated 
to  B.  BC  is  at  (267?2,  15'.'66)  in  2MASS. 

HIP  4878  has  faint  stars  B  and  C  at  2'.'9  and  32'.'2  distance,  re¬ 
spectively  (RAO  39).  The  optical  nature  of  C  can  be  established 
by  its  different  position  (145?3,  32'.'4)  and  its  color  in  2MASS. 
AB  was  remeasured  with  the  P200  one  year  later  and  found  at  a 
slightly  different  position,  so  it  is  likely  optical  as  well,  although 
the  field  is  not  crowded. 

HIP  5276  is  an  SB  with  P  =  49.5  d  (Gorynya  &  Tokovinin 
2014),  accompanied  by  the  newly  resolved  tertiary  component 
RAO  40  at  6'.'2.  The  physical  nature  of  the  tertiary  is  confirmed 
by  the  2MASS  position  and  by  observations  with  P200. 

HIP  5313  is  another  triple  system:  Aa,Ab  is  an  acceleration 
and  spectroscopic  binary  with  a  distant  physical  component  B 
(RAO  4 1 )  at  26'.'2,  confirmed  by  its  position  and  color  in  2MASS . 

HIP  6653  has  a  new  companion  at  1'.'8  (RAO  6).  The  star  is 
on  the  Keck  exoplanet  program,  its  radial-velocity  (RV)  trend 
(D.  Fischer,  private  communication,  2012)  could  be  caused  by 
the  companion.  The  companion  B  is  red,  the  field  has  low 
crowding,  so  the  physical  nature  of  AB  is  likely. 

HIP  7845  was  measured  as  a  test  target,  it  is  not  in  the 
FG-67  sample.  The  faint  companion  C  at  24"  moved  substan¬ 
tially  during  a  year,  hence  it  is  optical. 

HIP  11696  is  located  in  a  crowded  region  of  the  sky.  We 
targeted  the  component  B  and  found  a  faint  star  at  15'.'7  from  it 
(RAO  42BE),  considered  optical  (just  as  the  other  components 
C  and  D  listed  in  the  WDS)  because  of  the  high  crowding, 
non-hierarchical  configuration,  and  the  2MASS  position  of  BE. 

HIP  12062  has  an  optical  companion  RAO  43  at  14'.'8. 
The  pair  has  a  different  position  in  2MASS  and  the  region 
is  crowded. 

HIP  12067  has  a  new  5'.'7  companion  B  (RAO  7),  confirmed 
by  the  2MASS  position,  in  addition  to  the  inner  acceleration  and 
spectroscopic  binary  with  yet  undetermined  period. 


HIP  12189  and  HIP  12184  form  a  38"  physical  pair  STEA 

5  where  each  component  is,  in  turn,  a  close  binary:  A  is  SB 
(P  =  1.1  d),  B  is  resolved  here  into  a  0'.'53  pair  RAO  8, 
confirmed  next  year  with  Robo-AO  and  P200.  Moreover,  an 
exoplanet  with  P  =  335  d  around  B  was  announced.  We  discuss 
this  multiple  system  in  a  separate  paper. 

HIP  12685  has  a  companion  at  21'.'7  (RAO  44)  that  is  also 
present  in  2MASS  at  approximately  the  same  position.  However, 
the  PM  is  small,  the  field  is  crowded,  and  the  color  of  B  does 
not  match  a  dwarf  at  the  same  distance  as  A,  so  we  consider  B 
optical. 

HIP  12925  has  a  new  companion  at  1'.'9,  measured  several 
times  with  Robo-AO  and  confirmed  with  the  P200.  This  is  a 
young  multiple  system:  the  radial  velocity  (RV)  of  A  is  variable 
(Nordstrom  et  al.  2004),  there  is  a  CPM  companion  at  494" 
(Tokovinin  &  Lepine  2012),  and  HIP  12862  at  54'is  co-moving. 

HIP  13336  contains  the  visual  binary  A  281AB  and  a 
companion  RAO  45 AC  at  15",  considered  optical  because  it 
is  found  at  a  different  position  in  2MASS. 

HIP  15329  contains  the  AB  pair  (STE  53,  P  =  113  yr)  and  a 
CPM  companion  C  (7  =  13.24)  at  107"  found  by  (Tokovinin 

6  Lepine  2012).  Here  C  is  resolved  at  1'.'6  (RAO  10Ca,Cb), 
revealing  this  system  as  a  2-^2  quadruple.  We  believe  that  Ca,Cb 
is  physical  because  of  its  small  separation,  moderate  crowding, 
and  the  red  color  of  Cb  (A/  =  4.0  mag,  Az  =  2.8  mag). 

HIP  16329  Ba,Bb  =  RAO  llBa,Bb  (also  known  as  J  207)  is 
a  spectroscopic  binary  in  the  Hyades  cluster  with  a  period  of 
about  30.74  yr  (Griffin  2012).  Assuming  the  mass  of  the  K7V 
component  Ba  to  be  0.1  Mq,  the  spectroscopic  orbit  implies 
a  minimum  mass  of  0.3  Mq  for  Bb  and  a  semi-major  axis  of 
0'.'265.  The  pair  Ba,Bb  was  resolved  by  Robo-AO  three  times 
at  0'.'16  with  M  =  0.87'”,  Ar  =  0.97'”,  and  Az  =  0.52'”.  This 
corresponds  to  a  mass  for  Bb  of  ~0.6  Mq.  We  can’t  help  noting 
that  Bb  appears  redder  than  Ba,  and  that  its  lines  were  not 
detected  in  the  spectrum;  it  could  itself  by  a  close  binary  or  a  fast 
rotator.  This  pair  is  a  secondary  component  in  the  wide  binary 
STE  399AB  at  20"  separation,  also  measured  here  (despite  its 
WDS  designation,  AB  was  first  resolved  by  W.  Herschel  in 
1782,  not  by  W.  Struve).  The  primary  component  A  =  vB  3  is 
itself  a  close  binary  with  astrometric  acceleration  (Makarov  & 
Kaplan  2005),  although  Griffin  states  that  RV(A)  is  constant;  it 
is  not  resolved  by  Robo-AO. 

HIP  17022  has  an  astrometric  orbit  with  P  =  3.06 yr.  This 
binary  is  too  close  to  be  resolved  with  Robo-AO,  but  we  found 
another  faint  star  at  1'.'7  (RAO  47)  that  is  confirmed  as  physical 
with  the  P200. 

HIP  17217  is  a  triple  system  with  inner  binary  BU  1181  and 
a  new  component  at  4'.'8  (RAO  48AC)  confirmed  as  physical  by 
its  presence  in  the  2MASS  image. 

HIP  17336  has  a  new  companion  at  lOl'.A  that  is  not  found  in 
the  2MASS  point-source  catalog,  but  is  seen  at  270°  and  12" 
from  A  in  the  /-band  image.  The  difference  in  positions  matches 
the  reflex  PM  of  A,  hence  the  new  companion  is  optical. 

HIP  18218  is  the  binary  A  1293AB.  The  new  companion  C  at 
26"  is  optical,  as  evidenced  by  its  position  in  2MASS  and  high 
crowding. 

HIP  18719  is  an  SB  in  the  Hyades  (Griffin  2012)  that  is  not 
resolved  here.  The  faint  star  at  5'.'9  found  here  is  optical  because 
it  is  not  seen  in  the  2MASS  images. 

HIP  19389  has  a  CPM  companion  C  at  61"  (Tokovinin  & 
Lepine  2012),  targeted  here.  The  pair  CD  is  found  in  2MASS, 
which  shows  that  D  is  optical  to  C,  as  the  change  in  its  position 
simply  reflects  the  PM(C). 
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HIP  21099  is  a  member  of  the  Hyades.  The  faint  0'/9  pair  BC 
seen  at  19''  from  A  is  optical,  as  evidenced  by  its  position  in 
2MASS. 

HIP  21443  has  a  new  companion  B  at  5'.' 7,  revealed  as  physical 
by  its  color  and  position  in  2MASS  and  confirmed  with  the  P200. 
However,  the  PM(A)  is  small  and  the  region  is  crowded,  so  the 
status  of  the  new  companion  is  not  certain.  The  main  target  is 
also  an  SBl  with  a  period  of  2.06  d  (Gorynya  &  Tokovinin 
2014). 

HIP  23396  is  the  binary  HU  445  that  deviates  from  its  orbit. 
The  star  C  at  38"  is  classified  as  optical  by  comparing  with  its 
position  in  2MASS. 

HIP  24016  is  a  double-lined  SB  with  a  CPM  companion  B  = 
HIP  24005  at  69".  We  targeted  B  and  found  a  star  C  at  106?2, 
26'.'2  from  it  in  this  crowded  field.  The  star  is  at  103.9?  1  and 
27'.'32  in  2MASS,  hence  it  is  optical. 

HIP  25300  is  a  system  consisting  of  the  I'.'l  visual  binary 
STF  677AB  and  a  new  tertiary  component  C  at  6'.' 8,  confirmed 
as  physical  by  its  presence  in  the  2MASS  images.  Observations 
with  the  P200  also  confirmed  C  and  resolved  B  into  a  close  pair. 

HIP  26444  is  only  triple — a  single-lined  SB  coupled  to  the 
CPM  component  B  at  229"  distance  (Tokovinin  &  Lepine  2012). 
We  targeted  B  and  found  two  stars  around  it  at  "^4".  They  are  not 
present  in  the  2MASS  ^^-band  image  and  hence  are  considered 
optical. 

HIP  27067  and  HIP  27070  form  the  22"  binary  STF  775.  We 
measured  another  star  C  at  11"  from  A  in  this  crowded  field; 
AC  is  optical  judging  by  its  position  in  2MASS. 

HIP  27246  is  a  new  triple.  The  11.1  yr  spectroscopic  and 
astrometric  binary  was  recently  resolved  by  Horch  et  al.  (2012). 
The  new  component  C  at  \Q('.9  is  confirmed  as  physical  by  its 
position  and  color  in  2MASS. 

HIP  31267  is  another  triple  system  consisting  of  an  SB  and 
the  newly  found  tertiary  component  C,  confirmed  as  physical 
by  the  2MASS  image  and  follow-up  observations  with  the 
P200.  The  53"  pair  UC  1450AB  is  not  physical  because  the 
red  color  of  the  secondary  places  it  well  above  the  main 
sequence,  while  the  small  common  PM  is  likely  a  chance 
coincidence. 

HIP  33355  is  an  SB  with  a  new  companion  at  5'.'5  that  could 
be  physical.  The  ^^-band  image  in  2MASS  is  extended  toward 
a  position  angle  of  90°  and  the  companion  is  confirmed  with  the 
P200.  Yet,  the  small  PM  of  the  target  makes  it  difficult  to  make 
a  firm  conclusion  on  the  companion’s  status. 

HIP  34110  is  a  triple  system  with  a  5.1  d  SB  in  the  visual 
pair  A  619,  which  deviates  substantially  from  its  orbit.  The  new 
component  C  at  15'.'7  is  optical  as  shown  the  2MASS  data. 

HIP  35265  has  a  wide  CPM  companion  B  at  927"  (Tokovinin 
&  Lepine  2012),  whose  physical  nature  is  not  certain,  but  likely, 
given  the  substantial  PM.  We  targeted  B  and  measured  RAO 
58BC  at  17'.'6.  This  pair  has  a  different  position  in  2MASS, 
hence  is  optical. 

HIP  38018  has  an  optical  companion  C  at  118?7,  15'.'4,  with 
a  different  position  in  2MASS. 

HIP  40298  has  a  CPM  companion  B  at  244"  (LDS  2564), 
which  was  targeted.  The  pair  BC  found  here  is  optical. 

HIP  40479  has  an  optical  companion  at  31",  with  different 
position  in  2MASS. 

HIP  40918  has  a  CPM  companion  B  =  HIP  40882  at  258" 
which  is  resolved  here  into  the  2'.'8  pair  RAO  13BC  and 
confirmed  by  the  2MASS  image.  Considering  that  A  itself  is 
a  long-period  SB,  this  is  another  2-^2  quadruple  discovered  by 
this  survey. 


HIP  41319  has  a  CPM  companion  at  692",  likely  physical 
considering  its  large  PM.  It  was  targeted  here,  but  the  16'.'8  pair 
RAO  62BC  is  optical,  as  revealed  by  its  different  position  in 
2MASS. 

HIP  43426  is  a  visual  triple  with  two  companions  at  3'.'6 
and  49"  from  it.  The  distant  component  C  was  targeted.  We 
measured  RAO  63CD  at  18'.'9,  but  this  pair  is  optical. 

HIP  48273  is  a  double-lined  SB  with  P  =  3.955  d  (Griffin 
&  Suchkov  2003)  and  has  a  CPM  companion  B  at  1155" 
(Tokovinin  &  Lepine  2012),  considered  physical  because  of  the 
substantial  common  PM  and  matching  colors.  This  secondary 
was  resolved  here  at  O'.' 16  (RAO  90Ba,Bb),  converting  the 
system  into  a  quadruple.  The  pair  Ba,Bb  was  measured  at  SOAR 
on  2013.13  at  similar  position. 

HIP  49638  is  an  astrometric  binary  with  P  =  4.2  yr  and 
an  estimated  semi-major  axis  of  0'.'05.  We  detect  the  new 
companion  RAO  15  at  0'.'5,  converting  this  into  a  triple  system. 

HIP  69160  is  a  double-lined  SB  with  P  =  8.4  yr  (D.  Latham, 
private  communication,  2012)  and  expected  semi-major  axis  of 
0'.'09.  It  is  resolved  here  at  O'.' 12. 

HIP  69322  has  a  CPM  component  B  (V  =  13.92)  at  606" 
(Tokovinin  &  Lepine  2012),  which  is  resolved  twice  into  a  new 
0'.'58  pair  RAO  17BC. 

HIP  71843  has  a  CPM  component  B  at  216",  U  =  14.9.  We 
measured  a  pair  BC  at  17'.'3  twice.  Its  position  in  the  2MASS 
is  similar,  but  C  is  “bluer”  than  B,  hence  it  could  be  a  white 
dwarf,  although  more  likely  it  is  a  background  star,  considering 
the  small  PM  of  the  system. 

HIP  75676  is  a  40"  physical  binary  KU  108  known  since 
1893.  Both  its  components  are  SBs  (D.  Latham,  private  com¬ 
munication,  2012).  In  addition,  we  found  the  new  0'.'4  pair  BC, 
making  this  a  quintuple  system. 

HIP  79607  is  a  known  quintuple  system.  The  component  E 
(STF  2032AE,  E  =  HIP  79551  =  GJ  615.2C)  is  resolved  here  at 
^'4  (but  not  for  the  first  time:  Ea,Eb=YSC  152).  Another  star 
F  at  32"  from  E  (RAO  65)  appears  optical;  it  is  not  found  in 
2MASS. 

HIP  79629  has  a  CPM  component  B  at  107",  U  =  14.2.  It  is 
resolved  here  into  a  O'.' 14  pair  of  equal  stars,  converting  binary 
into  a  triple. 

HIP  81608  has  a  CPM  component  B  at  179".  The 
27"  pair  RAO  82BC  is  optical,  it  has  a  different  position 
in  2MASS. 

HIP  85042  has  a  CPM  companion  B  at  49"  according  to 
Raghavan  et  al.  (2010).  The  secondary  is  resolved  here  at  0'.'75 
(RAO  19BC),  while  A  itself  is  single. 

HIP  86642  is  a  double-lined  SB  (D.  Latham,  private  com¬ 
munication,  2012)  with  a  newly  discovered  tertiary  component 
at  2'.'2.  The  2MASS  image  does  not  resolve  this  pair,  but  it 
is  confirmed  with  the  P200,  where  the  inner  binary  was  also 
resolved. 

HIP  89207  is  a  2"  binary  AB  (A  2260).  Another  star  C  found 
here  at  17"  is  likely  optical,  as  the  field  is  very  crowded. 

HIP  91120  has  a  CPM  companion  B  at  621",  U  =  11.14.  It 
is  resolved  here  at  O'.' 14  (RAO  83Ba,Bb). 

HIP  94540  is  a  visual  triple  system  where  the  1"  pair  BU 
975BC  is  located  at  33'.'6  from  A.  However,  it  was  outside  the 
Robo-AO  field,  and  we  measured  only  the  known  optical  pair 
AD  at  15'.'9  and  the  new  companion  E  at  2'.'9  (RAO  84 AE).  The 
star  E  is  not  seen  in  the  ^^-band  2MASS  image  and  is  presumed 
to  be  optical,  until  proven  otherwise. 

HIP  94666  is  an  SB  for  which  we  find  a  distant  companion 
at  3'.' 6,  confirmed  in  the  2MASS  image. 
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HIP  94905  is  a  double-lined  SB  with  P  =  5.38  d.  The 
companion  RAO  67  at  6''9  has  the  same  position  in  2MASS. 
However,  the  field  is  extremely  crowded  and  the  PM(A)  is  small, 
so  the  new  companion  is  considered  to  be  likely  optical. 

HIP  95309  is  similar  to  the  previous  case:  we  find  a  new 
5'.'  1  companion  to  the  SB  in  a  very  crowded  field.  However,  the 
PM(A)  is  0'.'2yr“\  so  comparison  with  the  2MASS  position 
confirms  the  physical  nature  of  RAO  68. 

HIP  95769  is  a  spectroscopic  and  astrometric  binary  with 
P  =  2.26  yr.  The  companion  B  at  9^ '.9  is  optical;  it  has  a  different 
position  in  2MASS  and  the  field  is  crowded. 

HIP  96395  is  an  SB  converted  into  a  triple  by  establishing 
the  physical  nature  of  the  10''  companion  noted  previously  by 
Fuhrmann  (2004).  Its  position  in  2MASS  is  same  as  measured 
here,  its  color  matches  a  main-sequence  dwarf. 

HIP  97222  is  the  binary  STF  2576.  It  forms  a  quadruple  or 
quintuple  system  with  HIP  97295  at  792",  which  is  the  primary 
(Raghavan  et  al.  2010).  The  field  is  crowded  and  several  optical 
components  are  listed  in  the  WDS.  We  add  to  this  list  another 
optical  pair  RAO  69FK  at  16'.'4. 

HIP  99232  appears  to  be  triple.  It  was  resolved  at  0'.'065 
as  WSI  109  (not  confirmed  yet  by  a  second  measure),  while 
Robo-AO  finds  a  faint  companion  C  at  1'.'9,  reobserved  within  a 
year.  We  consider  RAO  2 1  AC  physical  (the  field  is  not  crowded), 
although  a  larger  time  base  is  needed  for  a  solid  confirmation 
because  PM(A)  is  small. 

HIP  99572  has  a  CPM  companion  C  at  1021"  (AC  =  TDT 
2085)  that  was  targeted  by  Robo-AO.  The  18"  pair  RAO  88CD 
is  optical  (D  is  not  found  in  2MASS). 

HIP  99965:  we  targeted  the  CPM  companion  F  (F  =  13.5, 
AF  =  GIC  155)  at  106"  and  measured  another  star  in  the  field; 
the  pair  RAO  70FL  is  optical,  its  position  in  2MASS  is  very 
different. 

HIP  101234  is  an  acceleration  binary,  resolved  here  for  the 
first  time  into  a  O'.' 17  pair  of  equal  stars  RAO  22,  and  confirmed 
with  the  P200. 

HIP  101430  appears  to  be  a  quadruple  system.  The  outer 
17"  pair  AB  =HJ  1535  is  composed  of  the  spectroscopic  and 
astrometric  binary  Aa,Ab  and  a  O'.' 17  pair  Ba,Bb  tentatively  seen 
in  the  pipeline  RoboAO  image,  but  not  confirmed  by  speckle 
processing.  This  subsystem  is  not  listed  in  Table  6  for  that 
reason,  but  it  was  resolved  with  the  P200  nevertheless.  The  12" 
pair  RAO  71AE  is  optical,  with  a  different  position  in  2MASS 
in  a  crowded  field. 

HIP  102040:  The  CPM  companion  C  at  125"  from  A  (AC  = 
LDS  1045)  was  targeted  and  resolved  at  0'.'27  (RAO  23Ca,Cb). 
The  two  other  wide  companions  D  and  E  are  optical;  they  are 
not  recovered  in  2MASS,  while  the  field  is  crowded. 

HIP  103455  is  an  acceleration  binary  resolved  here  at  0'.'6 
(RAO  24)  and  confirmed  later  with  the  P200.  Metchev  & 
Hillenbrand  (2009)  did  not  resolve  this  system,  although  its 
separation  implies  an  orbital  period  of  several  decades. 

HIP  103641  is  a  quintuple  system  that  consists  of  two  visual 
pairs  COU  2431Aa,Ab  (Aa  is  also  a  close  SB)  and  HDS 
2989Da,Db  (HIP  103052)  at  1132"  from  A.  We  targeted  D 
and  measured  the  close  pair,  as  well  as  an  optical  star  at 
19'.'7  from  it  (RAO  72DE),  which  has  a  different  position 
in  2MASS.  The  WDS  components  B  and  C  are  also  optical 
(crowded  field). 

HIP  104514  has  variable  RV  according  to  Nordstrom  et  al. 
(2004).  We  found  a  new  physical  companion  at  3'.'4  and 
measured  it  twice.  The  companion  is  confirmed  in  the  2MASS 
image  and  with  the  P200. 


HIP  108473  is  a  7.18  yr  SB  with  a  new  physical  companion 
RAO  73  at  12'.'4,  confirmed  by  its  position  and  color  in  2MASS. 

HIP  109361  is  an  acceleration  binary  with  variable  RV, 
resolved  here  at  0'.'37  (RAO  26). 

HIP  110291  is  astrometric  and  spectroscopic  binary  with 
P  =  2.11  yr.  The  two  wide  companions  B  and  C  measured  in 
this  crowded  field  are  optical,  as  evidenced  by  their  respective 
positions  in  2MASS. 

HIP  110574  is  acceleration  binary  resolved  here  at  0'.'09  and 
confirmed  with  the  P200  and  at  SOAR.  The  CPM  companion  B 
at  171"(y  =  15.88)  was  not  targeted,  being  too  faint. 

HIP  1 10626  is  also  an  acceleration  binary,  but  the  new  phys¬ 
ical  companion  at  4'.'4  discovered  by  Robo-AO  (and  confirmed 
in  the  2MASS  image)  is  too  distant  to  cause  the  acceleration. 
The  system  is  therefore  triple. 

HIP  111148  has  a  CPM  companion  B  (F  =  13.74)  at  62", 
AB  =  LEP  108.  We  targeted  B  and  found  two  faint  stars  C  and  D 
at  3'.'9  and  12'.' 3,  respectively.  Both  companions  are  considered 
optical  (C  is  not  seen  in  the  2MASS  image,  D  has  a  different 
position  of  268?3,  14'.'59,  and  the  field  is  crowded). 

HIP  112935  is  triple.  Its  CPM  companion  D  (F  =  13.3, 
AB  =  LDS  6388)  at  250"  is  resolved  here  at  ^'24  several 
times  at  fixed  position  during  one  year,  despite  the  fast  PM 
of  O'.' 5  yr“^ 

HIP  1 14456  has  a  CPM  companion  C  at  50"  (Raghavan  et  al. 
2010)  (the  WDS  pair  HI  1853AB  is  optical).  We  found  only  two 
optical  stars  around  C,  RAO  89CD  and  CE.  They  have  different 
positions  in  2MASS. 

HIP  115655  has  a  CPM  companion  B  =  HD  20748  at  185". 
We  measure  a  faint  CD  pair  at  13"  from  B.  The  color  of  C  in 
2MASS,  J  —  K  =  0.65,  suggests  that  it  is  optical,  although 
PM(A)  is  small  and  the  position  of  BC  in  2MASS  is  close  to  its 
measured  position. 

HIP  116906  has  a  CPM  companion  B  at  109'.'5  (LDS  5112, 
F  =  14.5),  which  is  resolved  here  into  a  0'.'5  pair  RAO  3  lBa,Bb 
(three  measurements  within  a  year).  The  main  star  hosts  an 
exoplanet  with  P  =  572  d. 

HIP  118213  is  an  acceleration  binary.  We  discover  a  tertiary 
companion  at  4'.' 8,  confirmed  by  its  repeated  measurement 
within  a  year  and  seen  in  the  2MASS  image.  The  new  pair 
RAO  76  was  also  measured  with  the  P200,  where  the  inner 
binary  was  resolved  as  well. 

HIP  118225  is  an  SB  with  P  =  25.4  d.  We  discovered  another 
component  at  5",  confirmed  by  its  repeated  measurement  and 
by  the  2MASS  image. 

Table  8  summarizes  the  position  measurements  for  the  Robo- 
AO  and  2MASS  data  used  to  determine  if  components  are  physi¬ 
cal  or  optical  pairs  in  many  of  the  newly  resolved  systems.  Some 
of  the  Robo-AO  pairs  were  not  resolvable  by  2MASS,  so  they  are 
not  included  in  the  table.  The  table  lists  the  binaries  by  HIP  num¬ 
ber,  and  the  components  that  are  being  measured,  followed  by 
the  measured  angle  (in  degrees)  and  separation  (in  arcseconds) 
for  the  Robo-AO  and  2MASS  data,  and  the  epoch  of  observation 
by  2MASS.  The  column  A  gives  the  total  displacement  between 
the  Robo-AO  and  2MASS  measurements.  The  next  column  ixAt 
gives  the  displacement  produced  by  proper  motion,  assuming 
At  =  12yr.  The  last  two  columns  are  the  proper  motion  of  the 
main  target. 

The  2MASS  measurements  at  "^5"  separations  are  likely  less 
accurate  than  its  astrometry  in  general  (~70  mas).  The  Robo-AO 
astrometry  of  wide  pairs  is  accurate  to  "^(X.'  1  (Law  et  al.  2014b). 
The  PM  of  companion  candidates  is  unknown  and  comparable 
to  the  PM  of  some  targets.  Considering  that  the  relative  PM  is 


19 


The  Astrophysical  Journal,  799:4  (21pp),  2015  January  20 


Riddle  et  al. 


Table  8 

Comparison  of  Positions  Measured  by  Robo-AO  with  2MASS 


HIPl 

Comp 

6*rao 

(deg) 

Prao 

(arcsec) 

2MASS 

Epoch 

6*2MASS 

(deg) 

P2MASS 

(arcsec) 

A 

(arcsec) 

/xAr 

(arcsec) 

< 

(masyr“^) 

(masyr“^) 

4016 

BC 

270.4 

17.00 

1998.97 

267.2 

15.66 

1.6 

1.2 

62 

-79 

4878 

AC 

147.0 

32.16 

2000.79 

145.3 

32.42 

1.0 

0.8 

66 

-23 

5276 

AB 

335.4 

6.17 

1999.97 

338.0 

5.32 

0.9 

1.6 

118 

-55 

5313 

AB 

244.7 

26.24 

1998.79 

245.0 

25.99 

0.3 

1.0 

-79 

-2 

7845 

AC 

175.3 

24.84 

2000.01 

174.6 

24.44 

0.5 

0.0 

-2 

2 

11696 

BE 

207.0 

15.68 

1999.95 

201.6 

16.68 

1.8 

1.5 

81 

-97 

12062 

AB 

207.8 

14.47 

2003.60 

199.4 

14.09 

4.2 

2.9 

230 

-83 

12067 

AB 

169.0 

5.71 

2000.89 

172.4 

5.25 

0.6 

1.5 

97 

-76 

12685 

AB 

211.1 

21.67 

1997.77 

212.2 

21.26 

0.6 

0.3 

-22 

4 

13336 

AC 

81.8 

15.20 

2000.75 

92.3 

20.79 

6.5 

6.8 

487 

-283 

18218 

AC 

241.6 

25.91 

1999.78 

241.4 

26.18 

0.3 

1.0 

-43 

-67 

19389 

CD 

294.3 

9.58 

1999.76 

284.2 

9.39 

1.7 

1.6 

-25 

-126 

21099 

AB 

310.3 

18.80 

1999.83 

310.9 

17.41 

1.4 

1.3 

102 

-41 

21443 

AB 

14.4 

5.74 

2000.00 

12.6 

5.32 

0.5 

0.2 

-5 

20 

23396 

AC 

242.8 

38.41 

1997.84 

242.1 

40.04 

1.7 

1.3 

-63 

-92 

24016 

BC 

106.2 

26.15 

1997.93 

103.9 

27.32 

1.6 

0.7 

58 

-5 

27067 

AC 

151.9 

11.10 

2000.86 

154.2 

12.07 

1.1 

1.1 

20 

-85 

27246 

AC 

68.5 

10.85 

2000.23 

67.9 

10.93 

0.1 

2.8 

142 

-187 

34110 

AC 

10.6 

15.71 

1998.94 

8.3 

15.23 

0.8 

0.7 

-50 

-20 

35265 

BC 

126.6 

17.64 

1999.98 

135.4 

18.43 

2.9 

2.2 

-58 

-173 

38018 

AC 

118.8 

15.48 

1999.08 

127.1 

15.66 

2.3 

2.1 

-93 

-151 

40298 

BC 

37.6 

13.07 

2000.17 

39.2 

15.53 

2.5 

2.2 

142 

113 

40479 

AB 

170.6 

31.25 

2000.07 

169.1 

32.09 

1.2 

0.9 

71 

-15 

41319 

BC 

310.0 

16.80 

1997.85 

298.4 

18.15 

3.8 

2.9 

-192 

-152 

43426 

CD 

357.6 

18.94 

1999.85 

350.9 

20.72 

2.9 

2.5 

-173 

113 

71843 

BC 

61.9 

17.27 

2000.20 

61.7 

16.38 

0.9 

0.7 

-58 

-18 

81608 

BC 

148.9 

27.16 

1999.27 

151.9 

28.53 

2.0 

1.2 

-28 

-98 

89207 

AC 

71.9 

17.24 

2000.54 

72.9 

16.51 

0.8 

0.5 

-1 

-39 

94905 

AB 

190.9 

6.95 

2000.32 

192.7 

6.68 

0.3 

0.4 

-29 

-10 

95309 

AB 

350.9 

5.13 

1997.45 

349.9 

4.80 

0.3 

2.2 

13 

-186 

95769 

AB 

348.2 

9.92 

1997.86 

341.8 

12.02 

2.4 

2.8 

49 

229 

96395 

AB 

250.7 

10.08 

1997.47 

249.3 

10.06 

0.2 

2.4 

0 

-200 

97222 

FK 

318.3 

16.40 

1998.33 

304.1 

13.05 

6.5 

5.3 

13 

-441 

99965 

EL 

264.2 

15.60 

2000.37 

257.3 

18.70 

3.7 

3.7 

-221 

-216 

101430 

AE 

224.7 

12.28 

1998.48 

224.4 

9.12 

3.2 

2.5 

154 

140 

103641 

DE 

117.6 

19.75 

2000.36 

108.4 

21.29 

3.6 

3.8 

232 

209 

108473 

AB 

89.3 

12.39 

2000.75 

89.5 

12.37 

0.1 

3.4 

-248 

-133 

110291 

AB 

122.9 

14.98 

1998.78 

131.5 

13.62 

2.5 

2.4 

-190 

-64 

110291 

AC 

232.1 

30.57 

1998.78 

233.1 

33.69 

3.2 

2.4 

-190 

-64 

111148 

BD 

264.4 

12.28 

1998.78 

268.3 

14.59 

2.5 

2.2 

-176 

34 

114456 

CD 

22.3 

20.59 

1999.76 

23.7 

16.84 

3.8 

3.6 

-86 

-286 

114456 

CE 

290.6 

15.33 

1999.76 

273.6 

15.74 

4.6 

3.6 

-86 

-286 

115655 

BC 

65.4 

13.20 

1998.80 

66.8 

13.81 

0.7 

0.8 

47 

-45 

Note.  The  measurement  for  HIP  12062  comes  from  Zacharias  et  al.  (2013),  all  others  are  2MASS  positions. 


not  known,  the  assumption  of  a  12  yr  time  base  is  good  enough 
for  estimating  the  expected  displacement. 
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